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Abstract: Using the gauge/gravity correspondence, we study the properties of 2-point 
correlation functions of finite-temperature strongly coupled gauge field theories, defined on 
a curved space of general spatial topology with a dual black hole description. We derive 
approximate asymptotic expressions for the correlation functions and their poles, supported 
by exact numerical calculations, and study their dependence on the dimension of spacetime 

^-j- and the spatial topology. The asymptotic structure of the correlation functions depends 

on the relation between the spatial curvature and the temperature, and is noticeable when 
they are of the same order. In the case of a hyperbolic topology, a specific temperature 

ON is identified for which exact analytical solutions exist for all types of perturbations. The 

asymptotic structure of the correlation functions poles is found to behave in a non-smooth 
manner when approaching this temperature. 
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1 Introduction 

The aim of this paper is to analyze the properties of 2-point correlation functions of finite- 
temperature strongly coupled (deconfined) gauge field theories defined on spacetimes with 
general spatial topology, and admit a dual black hole description via the AdS/CFT cor- 
respondence. Of specific interest is the dependence of these properties on the topology 
of the space on which the gauge theory is defined, and its temperature. We will see that 
the topology is indeed encoded in the properties of the correlation function and its poles, 
corresponding to quasinormal modes in the dual bulk spacetime: the sound of topology. 
In particular, the asymptotic "slope" of the poles encodes information about the spatial 
topology and the temperature, but is "universal" otherwise, i.e. it doesn't depend on the 
type of operators considered. 

One motivation for this work is the fact that important aspects of the strong coupling 
dynamics of gauge field theories can depend crucially on the topology of space-time. For 
instance, the confinement-deconfinement phase transition is a property of strongly coupled 
thermal conformal gauge theories defined on a space-time with spherical spatial topology. 
The theory is always in the deconfined phase in the planar and hyperbolic cases. Thus, 
one may hope that insights to the confinement-deconfinement phase transition mechanism 
can be gained by analyzing the strong coupling dynamics on different spatial topologies. 
This can be done by using the gravitational dual description. A second motivation is 
the recently studied correspondence between gravitational dynamics and strongly coupled 
condensed matter systems. One can naturally envision in and out of equilibrium condensed 
matter systems realized on different spatial topologies, where our study can be of much 
relevance. A third motivation, perhaps a more remote one, is the relevance of the interplay 
between topology and field theory dynamics to astrophysical/cosmological setups. 

We distinguish three classes of spatial topology : The spherical, the flat and the hy- 
perbolic (The relevant thermodynamic quantities and their dependence on the topology are 
given in Appendix B). Consider a gauge theory defined on a (d— l)-dimensional spacetime 
with the topology of R x ^5T 2 ' w here &d-2 ^ s a constant curvature (d — 2)-dimensional 
Riemannian manifold of arbitrary topology 2 . The two basic scales that characterize the 

In the spherical case, a Hawking-Page phase transition exists, which is interpreted in terms of the 
AdS/CFT correspondence as a transition between a higher temperature deconfined phase and a lower 
temperature confined phase (see [2], [3]). Phase transitions might exist also for the flat and hyperbolic 
cases, e.g. the transition between the black hole solution at a higher temperature and the Horowitz-Myers 
"AdS soliton" at a lower temperature (see [7]). 
2 We assume that d > 4. 
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theory are the temperature T and the scalar curvature of the ^5_2 manifold Rq T - In the 
analysis we will consider two different regimes: 

1. The regime of v T n — — )• 0: In this case the curvature of the f^_ 2 manifold may 
essentially be neglected, so we expect all physical quantities to behave as they do in 
the case of a flat (Rq T = 0) space. Therefore, in this limit one does not distinguish 
between different topologies (aside from the known changes of the Laplace operator 
spectra between different topologies). This limit includes the hydrodynamic limit 3 : 

L FT FT 

—?p > 0, where u denotes the frequency and L s is the spatial Laplace operator 

eigenvalue. 



2. The regime of vi y 1 = fixed: In this case it will be possible to distinguish between 
the different topologies. In particular we will study the limit of — Jr- = fixed and 
^ S> 1, where analytical asymptotic expressions will be derived for the quasinormal 
modes and the 2-point correlation functions. 

Note that the v ' ^ n > 1 regime does not exist in the black hole description in the case 
of spherical topology (since the temperature has a finite minimal value) - it corresponds to 
a different (confining) phase on the field theory side. Nor does it exist in the case of flat 
topology (for which Rq T = 0). While this regime does exist for the hyperbolic case, there 

is nothing separating it qualitatively from the v T n — = fixed regime in the context of this 
paper. For these reasons, we won't consider this regime separately here. 

The main new results of this work are (see Appendix A for a list of notations and 
definitions): 

• In the asymptotic limit, where — J*- ~ fixed and y>l, the following holds: 

— For all types of topologies and all the temperatures , the quasinormal frequencies 
(and therefore the poles of the field theory correlators) are of the form 

A n = A + nAA, (1.1) 

where AA depends on the topology and temperature. 

— In this limit we derive asymptotic expressions for the field theory 2-point re- 
tarded correlation functions of the operators dual to the different perturbation 
types. These are given in Equations 3.58 and 3.71 for the scalar perturba- 
tion modes, in Equations 3.74 and 3.77 for the longitudinal vector perturbation 
modes, and in Equations 3.81 and 3.84 for the transverse vector perturbation 
modes. 

• In the hyperbolic case there is a specific temperature T c , where the bulk "black hole" 
solution has no singularity, and is isometric to AdS. In this special case the following 
holds: 



3 The hydrodynamic condition > requires that ■ — > since the minimal eigenvalue 

is of the same order as ^J\R FT \. 
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— At T = T c we derive exact expressions for the field theory 2-point retarded 
correlation functions of the operators dual to the different perturbation types. 
These are given in Equations 3.23 and 3.28 for the scalar perturbation modes, 
in Equations 3.30 and 3.31 for the longitudinal vector perturbation modes and 
in Equations 3.33 and 3.34 for the transverse vector perturbation modes. 

— For T > T c , the asymptotic AA has a real component Re(AA) > 0. AsT-f T c , 
AX becomes imaginary according to Equation 2.71 for d = 4 (as previously 
calculated in [21]), and according to Equation 2.73 for d > 4. 

The paper is organized as follows. In Section 2, we will review some known results 
regarding quasinormal modes (formulated in our notation). We will calculate the QNM 
equations for massless scalar and vector perturbations, the exact solutions for the "spe- 
cial" case mentioned above (the hyperbolic case with the specific temperature T c ), the 
approximate quasinormal frequencies in the asymptotic limit (— fe- = fixed and y > 1) 
and the diffusion mode in the hydrodynamic limit. We will discuss the dependence of 
the asymptotic quasinormal modes on the topology and the temperature, in particular in 
the hyperbolic case around the temperature T c , and demonstrate it using some numerical 
results. In Section 3, we will use AdS/CFT dictionary to calculate the 2-point correlation 
functions of the dual field theory. We will obtain exact expressions for the hyperbolic case 
with T = T c and approximate expressions in the asymptotic limit for general topology and 
temperature. We will compare the asymptotic expressions with exact numerical solutions. 
Section 4 is devoted to a discussion. Details of calculations are outlined in the appendices. 

2 Quasi-Normal Modes 
2.1 QNM Definition 

Quasinormal modes of black holes or black branes are defined as the late-time oscillation 
modes of the black hole metric and the fields coupled to it, satisfying certain boundary 
conditions. Put differently, these are the eigenmodes of the linearized equations of motion 
over the black hole background. The boundary conditions are specified at the black hole 
horizon and at spatial infinity. The boundary condition at the horizon is chosen so that 
the solution corresponds to a wave ingoing into the horizon, while the boundary condition 
choice at spatial infinity depends on the asymptotic nature of the spacetime. 
Choosing the incoming wave solution gives the boundary condition: 



where z = corresponds to the BH horizon and z = 1 corresponds to spatial infinity. 
As for the boundary condition at spatial infinity, different boundary conditions have been 
investigated for different asymptotic spacetime geometries. In the context of AdS/CFT, 
the relevant boundary condition is the one that corresponds to the poles of the retarded 
correlator of the CFT operators dual to the investigated field. For the fields discussed here, 
this condition amounts to the Dirichlet boundary condition (see [13], [17] and [22]): 



~ z c 




= o . 
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2.2 QNM Equations 

We consider the d-dimensional black hole metric given by a metric of the form: 

1 



rls 2 
as bulk 



-f(r) dt A 



f(r) 



dr 2 + r 2 dO 2 



d-2 i 



(2.3) 



where dTt 2 d _ 2 = £) f j{9n)ij dx l dx^ is the inner metric of a constant curvature, (d-2)- 
dimensions Riemannian manifold with of same topology as the ^-2 manifold. 

From the vacuum Einstein equation with a cosmological constant one can deduce the 
form of f(r) to be 



Sir) = k + 



R 2 



„d—3 



(2.4) 



where R is related to the cosmological constant by R 2 



(rf-2)(rf-l) 



2A 



and k is related to 

the scalar curvature Rq of the manifold VLd-2 by (see, for example, [18]) k = (d-2)(d-3) • 
In terms of these coordinates, the horizon of the black hole is located at r = r + where 
S(r+) = 0, while the boundary is located at r — > oo. 

After making the transformation to the coordinate z = 1 — — , the metric takes the 



form: 



ds 2 



£ v , dr + 



p 2 5(z)(i-z) 2 



dz" 



dO 



d-2 ; 



where p = ^ , K 



R 



4- and 



1 + A-fl 



(1 + K)(1 



\d-l 



(2.5) 



(2.6) 



In terms of the new coordinate z, the horizon is located at z = while the boundary is 
located at z = 1. 

Let us assume that the spectrum of the scalar and vector Laplace operator Aq 4 defined 
on the manifold Qd—2 1S given by 



A n H L 2(x) = L 2 s H L 2(x) 
A n A L 2(x) = L 2 v A L 2{x) 



(2.7) 
(2.8) 



respectively. In order for a QNM corresponding to the frequency oj and the above Laplace 
operator eigenvalue to exist, the corresponding boundary conditions problem should have a 



non-trivial solution. Introducing the dimensionless parameters A = — — . 



UJr + n = he- n = Lu. 
-j?-, Hs — p , <lv — p 

and C = g'(0), we find the following equations (see Appendix C for a detailed derivation): 



For a (massless, minimally coupled) scalar QNM, the equation is 



n 



\d-2 



d z 



9{z) 



(1 



\d-2 



d z ifj 





[ A2 2 1 


+ 


.9{z) 



rp = 



along with the boundary conditions: 

^=0 



~ Z c 



z=l 







(2.9) 



(2.10) 



4 In the context of this work we shall define the Laplace operator as: A = Sd where d and S are the 
exterior derivative and codifferential operators respectively. 
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For a "longitudinal" vector QNM, the equation is 

1 



d z 



~g(z)(l-z) d - 4 d. 



.(1 -z) d -* 
along with the boundary conditions: 



4' 





\ A 2 2 1 


+ 






_9{z) 



r/> = 0, 



(2.11) 



z=0 



_ i\ 

~ Z c 



(i - z) d ~ A d z 

or equivalently the equation 



(1-z) 



d-4 



2 = 1 



, (2.12) 



2=1 



:i - z ) d - 4 d 



9{z) 



(l-z) 



d-4 



+ 



:9{z)d z g{z)d z %l) + 



A 2 - q 2 s g{z) 
along with the boundary conditions: 

• For a "transverse" vector QNM, the equation is 



^ = 0, (2.13) 



(2.14) 



fl - z) d - 4 d. 



9{z) 



(l-z) '- 4 
along with the boundary conditions: 

Vt=o 



d z i(j 





f A2 2 1 


+ 





rp = 



(2.15) 



(2.16) 



It can be easily seen from the above equations that if ip(z) is as a solution for any of the 
equations (satisfying the corresponding boundary conditions) with frequency A, then ip*(z) 
is a solution for the same equation (with the same boundary conditions) with frequency 
—A*. This means that the QNM frequencies are always symmetric with respect to the 
imaginary axis. 

A detailed derivation of the QNM Equations can be found in Appendix C. 



2.3 Exact solutions for the K 



1 case 



In the case where the £ld-2 manifold is hyperbolic, we consider the special state corre- 
sponding to the temperature: 



T c = 



1 



I rft\ 



2-kR 2ir\ (d-2)(d-3) 



(2.17) 



for which K = — 1, and the horizon radius is given by pc = y \k\. In this case, the QNM 
equations (Equations 2.9, 2.11 and 2.15) can be solved analytically. This has been done 
in [19], [20] and [12]. Here we summarize these solutions in our notation. 
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Making the transformation: 

w = g(z) = l-(l-z) 2 (2.18) 
rp = u> 7 (l -w) s <f> (2.19) 

(where 7 and 5 are defined in Appendix E for each perturbation type), one obtains the hy- 
pergeometric equation. The solution to the equation satisfying the incoming-wave bound- 
ary condition at the horizon is 2 -Pi (o-,b;c;w), where a, b and c are given in Appendix E 
for each case. 

Applying the QNM boundary condition at the AdS boundary we get for the scalar and 
transverse vector cases: 

2 F l(o , 6;C ;l)=^I<9^=0, (2.20) 
1 (c — a)L (c — 0) 

and for the longitudinal vector case (and using known formula fot the hypergeometric 
function, taking into account c — a — b = ^y^) we get: 



(1-w) 2 d u 



Fi(a,b;c;w) 

2 2 

(c — a)(c — b) 
c 2 

(c - a)(c - 6) r(c + l)r(c - a - b + 1) 



«;=i 

F^a^c+l;!) 



r(c-o + l)r(c-6 + l) 



. (2.21) 



The solution for which is given bye — a = —n or c — b = —n where n £ Z and n > 0, so 
that the quasi-normal frequencies are: 

• For the massless scalar case: 



Xn = -^{d+l)±Jq 2 -^{d-3) 2 -2ni . (2.22) 



For the longitudinal vector case: 



X n = - Z -(d-3)±Jq 2 - J(d-3) 2 -2m . (2.23) 



For the transverse vector case: 



A„ = -~(d - 1) ± \jq 2 - \{d - 5) 2 - 2m . (2.24) 
2.4 QNM Asymptotics (for K > -1) 

For the case K > — 1, one may find an approximate analytical expression for the n-th QNM 
frequency for n — > 00 - that is, in the asymptotic limit of — |^ = fixed and ^ S> 1. This 
has been done in [10] using the "monodromy" method. Here we summarize this calculation 
using our notations. The general stages of the calculation are as follows: 
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1. We bring each of the QNM equations to a Schrodinger form, with some effective 
potentials, using appropriate coordinate and function transformations. 

2. We find the poles of these effective potentials on the complex plane, and write ap- 
proximate forms for the equations in the vicinity of the poles. 

3. We solve the approximate equations around each pole using Bessel functions. 

4. Using the asymptotic limit of |A| 3> 1, we replace the Bessel functions with their 
asymptotic forms. 

5. We match the solutions on each region of the complex plane, and apply the appro- 
priate boundary conditions, to get equations for values of A that allow for non-trivial 
solutions. 

6. We solve the equations analytically and obtain an asymptotic expression of the form: 

A n = A + nAA . (2.25) 

2.4.1 Calculation of Asymptotic QNMs 

Starting from the QNM Equations, we perform the following transformation: 

i/>(z) = (l-z)%<f>(z) (2.26) 

^ d-1 

dz = y~] 7fc ln(z - z k ) + z , (2.27) 
9{z) ^ 



k=l 



where: 



[ d — 2 for scalar, 

a = I , (2.28) 
[ d — 4 for vector 

Zk (k = 1, . . . , d — 1) are the zeros of the polynomial g(z) and 5 

z — Zk 1 , 

lk = lim — * = — - 2.29 
z-+z k g(z) g'{z k ) 

d-1 

^o = -^7fcln(l-z fc ) (2.30) 

k=l 

(so that z(z = 1) = 0). After this transformation the equations take the Schrodinger form: 

dj<t>+ [X 2 -V(z)} ( f> = , (2.31) 

where the effective potential V{z) is given for each case in Appendix G. 

Next we turn to the calculation of the asymptotic QNMs using the monodromy method. 
We first note that the asymptotic QNM frequencies satisfy (see [9] and [10]) 

Im(A5 ) = . (2.32) 



3 The branch cuts are chosen here so that < arg(z) < 2-k. 
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Looking at the anti-Stokes lines, defined by 



Im(AS) = 



(2.33) 



we choose two of them to be the contour for the calculation 6 . We then proceed by devel- 
oping the equation, and finding its solutions, around points along the chosen contour (the 
boundary and z — > oo ) and "matching" the solutions. The chosen contour is shown in 
Figure l 7 . 

We first look at the equation at z — > oo or z — > zq. As shown in Appendix G, the 
equation at S — > zq is approximately 



A 



Joo 



1 



. 



4(S-S ) 2 . 

The solution to this equation can be written as: 

« A + P + {~z - S ) + A^P^(z - zq) . 
In the non-scalar case where joo 7^ (and more generally joo isn't an even integer) 
P±(z - S ) = \/2tt\{z - So) J, joo (A(S - S )) , 



(2.34) 



(2.35) 



(2.36) 



where J, joo are the Bessel functions. The fact that asymptotically 8 |A| S> 1 allows us to 

=t 2 

replace the Bessel functions with their corresponding asymptotic forms. For direction (1) 
(as appears in Figure 1), we have A(S — So) > 0, so that 



A e i(-\za-a + ) _|_ ^_ e i(-A2o-a-) 



where 



+ 



7T 



-iXz 



(2.37) 
(2.38) 



For direction (2) we have A(S — So) < 0, so that 



A i(-\z +3a + ) _|_ £_ e i(-\z +3a-) 



iXz 



A +e i (^*o+a+) _|_ 7 4_ e «(Aio+a-) 



-iXz 



(2.39) 



Next we look at the equation at the boundary. As shown in Appendix G, the equation 
at S —7- is approximately 



A 



4(5 - z ) 5 



. 



(2.40) 



6 The choice of branch cuts so that < arg(z) < 2tt makes sure that the contours don't pass through a 
branch cut. Choosing < arg(z) < 2tt instead would give the complex conjugates of these contours. 
7 The assumed shape of the anti-Stokes lines is only true in the case of K > — 1. 

8 Another assumption is that |A| ^> |<j|, so that the effect of the effective potential can be neglected in 
the regions between the the "special" points (2 — > 0, z — ¥ 1, Z — ¥ 00) 
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Im(z) 
▲ 



z 



(2) 




(a) z plane contour 



lm(z) 

z 




(b) z plane contour 

Figure 1. Chosen contour for the monodromy calculation: (a) In the z plane and (b) In the z 
plane. 

The solution is therefore: 

4> « B + P + (z) + B^P^(z) = B + \/2ir\zJj 1 (\S) + y/2ir\zJ_ a (A5) . (2.41) 

2 2 

We now apply the boundary conditions at the boundary (z = 1, 5 = 0) to this solution. 
In all cases (the scalar, transverse vector and longitudinal vector) we end up with the 
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condition 



B_ = 



(2.42) 



Next we again assume |A| 3> 1 so that the Bessel function can be replaced with its asymp- 
totic form. For direction (2), we have \z > 0, so that 



where 



2B + cos(Az - P+) 



ft 



B+e 



-04 



e iXz + 



B+e 



i/S+ 



-iXz 



7T 



(i + ii) • 



(2.43) 



(2.44) 



Finally from the boundary conditions at the horizon (z = 0, z — >■ — oo) we have (for 
line (1) as it appears in Figure 1) 

<(> ~ e" iAi . (2.45) 

Now we "match" the solutions on lines (1) and (2). In the non-scalar case, in which 
ioo 7^ so that a+ 7^ a_. In this case we have for line (1) (from Equations 2.37 and 2.45) 



^ e i(-Ai -a+) _|_ ^_ e «(-Ai -a-) _ Q 

For line (2) we have (from Equations 2.39 and 2.43) 

A +e i(-\z +3a+) _|_ ^_ e i(-Af +3a_) _ Q +e ~^+ 
A+e i(Xz +a+) + A _ e i(Xz +a-) = B+ jP+ _ 

The condition for having a non-trivial solution for this system of equations is: 

e i3a+ s i n ( a+ ) _ e «3a_ s i n ( a _) 



tan(A^o — /3+) 
Finally, we obtain the result 9 : 



e i3a + cos ( a , + ) _ e «3«_ C os(a_ 

A„ = Ao + nAA , 



where 



/3_i_ + arctan 



e i3a + sin(a + )-e i3a - sin(q_) 
e i3a + cos(a + )-e' 3a - cos(a_) 



Z0 

In the scalar case, in which = 0, we have: 



AA = — 

zo 



P + (z - Zq) = v / 2vrA(z - z )Jo (A(z - z )) 
P_ {z - zq) = y/2ir\(z - z )Y (X(z - z )) 



(2.46) 



(2.47) 
(2.48) 



(2.49) 
(2.50) 
(2.51) 



(2.52) 
(2.53) 



In a calculation similar to the non-scalar case, we again obtain the asymptotic form in 
Equation 2.50, where 



A, 



P+ + f - § ln(2) 







AA 



7T 



Z() 



(2.54) 



9 Obviously, —A* would also be a solution. Choosing the branch cuts so that < arg(z) < 2n would give 
this solution instead. 



- 10 - 



2.4.2 Asymptotic Dependence on Topology and Temperature 

From the results of the calculation in this section we see that the asymptotic parameters of 
the QNM frequencies (namely Ao and AA) of each perturbation type depend on the £ld-2 
manifold topology and the temperature through the parameter zq, which depends only on 
d and K (where K contains the dependency on both the topology and the temperature). 
In fact, the asymptotic slope of the frequencies (AA) depends only on the spatial topology 
and the temperature - it is the same for all types of perturbations (even for types that 
aren't discussed here, such as tensorial perturbations): 

AA = — . (2.55) 

The dependence of z$ on K is plotted in Figure 2 for d = 4 and in Figure 3 for d = 5. In 
the case of spherical (K > 0) and flat (K = 0) topologies, the parameter zq, and therefore 
the slope AA, always has an imaginary part and a real part (for all temperature values). 
However in the hyperbolic case, at the temperature T c (as defined in Section 1) where 
K = — 1, the slope becomes completely imaginary 10 (as seen in the exact results given in 
Subsection 2.3 for this case). Moreover, it can be clearly seen in the Figures 2 and 3 that 
zq(K) (and therefore AA(T)) is not smooth at this temperature. 



!o(K) d=4 




Figure 2. The dependence of the parameter zo on the parameter K on the zq complex plane, in 
d = 4 dimensions. 

In the following we study two specific cases: the case of K = 0, and the case of K — > —1 
{k<0,T^T c ). 

The case of K = Let us first review the case of K = (where the £ld-2 manifold is 
"flat"), described for example in [10] as the limiting case of large spherical black holes. In 
this case, the roots of g(z) are analytically known and we can get an accurate expression 
for the asymptotic QNM frequencies (and most importantly the gap AA). 

10 In dimensions d = 4 and d — 5, Re So > below the temperature T c , so that the asymptotic slope of 
the QNM frequencies remains imaginary below this temperature as well. 
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Zo(K) d=5 





Figure 3. The dependence of the parameter zq on the parameter K on the zq complex plane, in 
d = 4 dimensions. 



Assuming that K = 0, we get for g(z): 

~g{z) = 1 - (1 - zf- 1 . (2.56) 

Its roots are given by: 

Zfc = l-e i2 ^H k = l,...,{d-l) . (2.57) 

After calculating the residues and using the definition of So, we obtain the results: 

vr 7t 1 

z = -- — -e t r (2.58) 



d-1 



sin 



d-1 



AA = — = -(d - l)e'«»-i sin ( 

z \d-l 



(2.59) 



For example, for d = 4 we get AA = — — if, and for d = 5 we get AA = —2(1 + i). 

The case of -PC — > — 1 Here we consider the case of K — > —1 (T — > T c ) from above 
(meaning K > —1). Define: 

K = -l + e e > , (2.60) 



so that 



5(Z) = 1-(1-6)(1-Z) 2 -6(1-Z) 



d-1 



(2.61) 
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Notice also that from the discussion in Section 1 and Appendix B we can relate e to the 
temperature difference by: 

9 AT 

(2.62) 



d-2 T r 



We first find an approximation for the roots of g(z). For K = —1 (e = 0) we have: 



g(z) = l-(l-zf 







z 1 = 0, z 2 = 2 . 



(2.63) 



In this case there are only two roots. When K > — 1, there are d — 1 roots. Two of them 
will be close to the two roots of the K = — 1 case, while the others go to oo as K — )■ — 1. 
Developing the roots and residues around z = and z = 2 to first order in e, we get: 



zi = 
7i 



f 



1 d-3 



2 + (d-3)e 2 

l + (-l) d 
z 2 ^2+ \ ' e 



72 



1 (-l) d (d-2) + l 



-e . 



Approximating the other roots we get: 

1 — Zk ~ e~^3e J ^ x « J 

-f <i-3e \d-3^ d-3 J 1 



7fc 



d-3 



(2.64) 
(2.65) 

(2.66) 
(2.67) 

(2.68) 
(2.69) 



where j = d — k — 1 = 0, . . . ,d — 4. 

We put these values into the definition for z$, and after some calculations we obtain 
the following results: 
For the case d = 4: 





v'tt 


+ e In e 


Z PS 


T 




7T 


4 , 


AA = 




ps -2i + -elne , 




^0 


7T 



And for the case d > 4: 





?7T 


7T 


1 i 




Y 


ri- 




AA = 


7T 




4 1 i 




ps -2% 


gd-3 




5o 




d " 3 « in (^3) 



For example, for d = 5 we get 



AA 



-2% - 2Je . 



(2.70) 
(2.71) 

(2.72) 
(2.73) 

(2.74) 
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2.5 Hydrodynamic Approximation 



, /\R FT \ 

Taking the hydrodynamic limit of X,q — > (y, ^ > 0, v T n > is required) in the 

appropriate QNM equations, one may calculate the hydrodynamic constants of the dual 
gauge theory (see [14], [16]). Since the longitudinal vector mode is coupled to charge 
density fluctuations in the gauge theory, its QNM spectrum in the hydrodynamic limit has 
to contain a mode corresponding to the charge diffusion mode in the gauge theory. Here 
we review the derivation of the diffusion mode from the QNM Equation in our notation. 
Starting from Equation 2.11 for a longitudinal vector QNM, let us define a new function 

i X 

if)' = z~~cip, and get a new "shifted" equation: 



0, 



iX 
~Cz 



cZ-4 



C~z 



1 



(1-z 



id— 4 



+ 



V> = 0, (2.75) 



with boundary conditions 

^| 2=0 ~ const. (1 - z) d ~ 4 ( 



Cz 



1 - z) d - 4 



. 



(2.76) 



2=1 



For the hydrodynamic approximation we shall assume A ~ q 2 <C 1, and develop tp in 
orders of A : 

V; = V (0) + + ^ (2) + • • • • (2-77) 

We put this into Equation 2.75, equate each order to and solve the equation with the 
corresponding boundary conditions in terms of the higher order solutions. 
For the 0-order we get the general solution: 

1 



V> (0) ~ C 



C 



In z + const. 



+ D (l-z 



,cZ-4 



since g(z) 



so that 



2=0 



Cz. From the boundary condition at z = we get 

C = , 



^(0) 



Dn(l 



\d-4 



(2.78) 

(2.79) 
(2.80) 



while the boundary conditions at z = 1 are automatically fulfilled. 
Next we turn to the A 1 order, for which we get the solution: 



(i) 



iXDo (l- z )*-*]nz-(l-z) d -* q2Do 



C 



+ d(l-z] 



d-3 

d-A 



1 



i 



dz 



g(z)(l-z) d -* 
Applying the boundary condition at z = we get 

i\D 



dz + D^l-z 



,d— 4 



2=0 



C 



— — £^ m z + — In z + i^i ~ const. , 



(2.81) 



(2.82) 
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from which we get the condition 



,2 



iX - f-^ ) D + d = . (2.83) 



From the boundary condition at z = 1 we have 



2 = 1 U 



2=1 



= ^%,4 + Cl - ^%, 4 = d = . (2.84) 
Putting this into Equation 2.83 we get to the conclusion 

»— sV- (2 ' 85) 

which is a (normalized) diffusion relation, with the diffusion constant equal to ^tz^, ov, in 
terms of u; and L^ 7 : 

W «-iZ)(Lf T ) 2 £> = ^t4- (2-86) 

This is true for any dimension d > 4 and any f2d_2 topology and temperature for whom 
the hydrodynamic condition can be fulfilled. As expected, the leading term in this approx- 
imation doesn't depend on the topology of 

2.6 Numerical Calculation of QNMs 

In the following we present some results of exact numerical calculations of the QNM spectra 
in several cases in order to illustrate the dependence on topology and temperature. The 
spectra have been calculated using the method outlined in Appendix I. 

Figures 4 and 5 contain the exact QNM spectra for the cases of spherical {K = 1), 
flat (K = 0) and hyperbolic {K = —0.5) topologies, for the scalar and longitudinal vector 
cases. The dependence of the asymptotic slope of the spectrum on the parameter K can 
be easily seen in these figures (as well as its independence of other parameters such as the 
spatial mode or perturbation type). The (pure imaginary) hydrodynamic mode can also 
be seen in the spectra of the longitudinal vector perturbations. 

The numerical results for —1.1 < K < —0.9, as illustrated in Figures 6 and 7, demon- 
strate the phenomenon discussed in Sections 2.3 and 2.4 for an hyperbolic £ld-2 manifold: 
At T = T c (K = — 1) the asymptotic QNM frequency gap (AA) becomes imaginary, and 
for d = 4, 5 it remains imaginary for T < T c (K < — 1). 
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Scalar d =4 q=0. 




Figure 4. Numerically calculated QNM fre- 
quencies for the case of scalar perturbations with 
d = 4,K = 1, 0, -0.5 and q s = 0. 



LonBVecliir d=5 q=1.41 




Figure 5. Numerically calculated QNM fre- 
quencies for the case of longitudinal vector per- 
turbations with d — 5,K — 1,0,-0.5 and q s = 
V2. 



Scalar .1=4 q=l. 




Figure 6. Numerically calculated QNM fre- 
quencies for the case of scalar perturbations with 
d = 4,-1.05 < K < -0.9 and q s = 1. 



l.im^V cclor d-S q — 1.5 




Figure 7. Numerically calculated QNM fre- 
quencies for the case of longitudinal vector per- 
turbations with d = 5,-1.1 < K < —0.9 and 
q s = 1.5. 
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3 CFT Correlators 



3.1 CFT Correlators From Holographic Principle 

Via the AdS/CFT corresondence, a black hole (or black brane) background in the bulk 
spacetime, such as the ones discussed here, corresponds to a dual CFT at finite temperature 
in the deconfined phase. The correspondence allows one to calculate the 2-point correlation 
function of CFT operators by calculating the on-shell action of the bulk field dual to that 
operator. 

The general prescription for calculating the dual CFT correlators from the AdS / CFT 
correspondence is discussed in [1], [2], [4], and [5]. In the case of Minkowski spacetime, the 
calculation of the Minkowski correlators from AdS/CFT requires some additional subtleties, 
as explained in [13]. In particular, the retarded green function defined as: 

G R (t, x; t', x') = -id(t - t') < [0(t, x), 0(t', x')] > (3.1) 

can be calculated by choosing the incoming-wave boundary condition at the black hole 
horizon for the solution of the classical EOM in the bulk spacetime (see Subsection 2.1). 

One consequence of the above prescription is the fact that the QNM frequency spec- 
trum of the bulk field (as defined in Subsection 2.1) comprises the poles of the retarded 
correlator of the dual CFT operator. Therefore all of the results of Section 2 serve to teach 
about the poles of the retarded correlators of the theories dual to the discussed black holes. 
The purpose of this section is the expansion of those results from knowledge of the poles 
to expressions for the correlators themselves. 

3.2 General Formulae For Correlators 

The background spacetime metric on which the field theory dual to the bulk black hole is 
defined will be given by: 

ds 2 FT = -dt 2 + dO| Td _ 2 , (3.2) 

where d£l 2 FT d2 = J2i j{9n T )ij dx l FT dx FT is the inner metric of the O^j^ manifold. This 
metric is conformally related to the bulk background metric at its boundary, with the 
conformal factor chosen to cancel the divergence in the bulk metric. The relation between 
the metrics is then given by (see [23]) 

C = }™(^9™ k ) , (3-3) 

or: 

ds| T = - dt 2 + R 2 dn 2 d _ 2 . (3.4) 

Comparing this with Equation (3.2) we get the relation between the field theory manifold 
and the bulk manifold dQ FTd _ 2 = R 2 dQ d _ 2 , and therefore the relation between the 

corresponding curvatures - Rq T = and k F x = (d^2)(d-3) = W> anc ^ ^ e re l a ti° n 
between the corresponding Laplace operator eigenvalues Lf T = 

Using the AdS/CFT prescription for a massless scalar bulk field and a gauge vector 
bulk field, we may find expressions for the retarded correlation functions of the dual field 
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theory operators. This is done by integrating by parts the expressions for the action of 
these fields, using the EOMs they satisfy to get an expression that depends only on the 
boundary values of these fields and taking the derivative of the expression with respect to 
the the boundary values. We find the following expressions: 



For a massless scalar field, the correlation function is given by: 

„rf-i 



G r (u,b) 



-2C„ 



1 



R d il-z) 



(3.5) 



Z->1 



where ^) u ,a satisfies the scalar QNM equation (Equation 2.9) along with the boundary 
conditions: 



~ z c 



2=0 



1 



z-H 



(3.6) 



C s is the appropriate normalization constant for the bulk scalar field (see Appendix 
C). 

For the longitudinal component of a vector field, the correlation function is given by: 

„d-3 



Gg(w,s) = 4C v r + d _ 2 (1 _^-4 ^,s 



(3.7) 



G^(cj,s)=i-G R (cj,s) 



Gjj,( W ,s) = -^G£(a;,s) 



(3.8) 
(3.9) 



where ^a, jS satisfies the longitudinal vector QNM equation (Equation 2.11) along with 
the boundary conditions: 



2=0 



Z ° Xw,b| z -*i = -2»(«)(1 " ^ ^2 



(1 - z) d -* 



1 . 



Z->1 



(3.10) 

G„ is the appropriate normalization constant for the bulk vector field (see Appendix 
C). 

For the transverse component of a vector field, the correlation function is given by: 



G£l(w,v) = -4G, 



(3.11) 



2->l 



where ^ w ,v satisfies the transverse vector QNM equation (Equation 2.15) along with 
the boundary conditions: 



2=0 



1 



z-»l 



(3.12) 



Several notes about taking the limit z — > 1: 
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1. The expression d z ^ tS (z) in the limit z — > 1 is calculated by normalizing ifj at z = 1— e, 
differentiating it and then taking e — > 1, so that: 



(3.13) 

Z->1 



2. When taking the limit z — > 1 , one must drop the contact terms - the terms that are 
polynomials in uj and L s and diverge as z — )■ 1 (These terms are removed by process 
of renormalization) . This can be done by developing d z ijj Uta (z) in orders of 1 — z, and 
then dropping all terms up to the lowest non-contact-term order. 

Note also that the spacetime coordinates correlators for the vector perturbations are 
related to the above longitudinal and transverse components by the following relations: 

G«(t, X ;t',x') = [ p^G& U ,B)H™(z)nFWe-*>V-* (3.14) 

s 

Gfax-M = [ ^ £ Sg (cj, s)H™ (x) ^[l^ (3.15) 
Gfj{t,x;t', x ')= [ ^£Gf„(c 

J s 

A detailed derivation of these expressions can be found in Appendix D. 
3.3 Exact Correlators for the K = — 1 Case 

As in Subsection 2.3, in the K = — 1 case (where the £ld-2 manifold is hyperbolic and 
T = T c ), analytical expressions for the CFT correlators may be found. This section is 
dedicated to the calculation of the correlators in this case. 

3.3.1 Scalar Correlator 

Continuing from the transformation defined in Appendix E, the expression in Equation 3.5 
can be written in terms of <b and w. 



i U- • O J ^ prp ^ prp e 



C n U.s) If/ 



R d (i- w ) d -^ <f> Ufi 



(3.17) 

W— ¥l 



(where a contact term has been dropped). 

As explained in Appendix E, the solution to the EOM with an incoming- wave bound- 
ary condition at the horizon is (p u ^ s = 2-F1 (a, b; c; w), where a, b and c are given by Equa- 
tions E.13 and E.14. Define: 

d — 1 , 
A = c-a-b= — ^— . (3.18) 

There are now two possible cases: 
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1. d is even, so that A is non-integer. In this case, the hypergeometric function satisfies 
the following connection formula: 



</>u, a ( w ) = 2-F1 (a,b;c;w) 



r(c)r(A) 



T(c-a)T(c-b) 

:(c)r(-A 
r(o)r(b) 



2 F l (a,b;l-A;l-w) 



+ r L A) (1 ~ 2*1 (c- a, c - b; A + 1; 1 - w) 



Defining the coefficients: 



we get: 



A(lj,s) 
B(uj,s) 



A g(^s) e A-l 



r( c )r(A) 



r(c-a)r(c-6) 

r(c)r(-A) 



«j=1— e 



A(w,s) 



r(a)r(6) 



A r(-A)r( c - Q )r(c-6)) cA _ 1 



r(A) r(a)r(6) 



(3.19) 

(3.20) 
(3.21) 

(3.22) 



Putting this into Equation 3.17 we get: 



r^" 1 r(-A)r(a + A)r(6 + A) 



G r {oj,s) = 2(d-l)a 



R d r(A) T(a) 



r(6) 



(3.23) 



2. d is odd, so that A is an integer. In this case, the hypergeometric function satisfies 11 : 
<Pui, s ( w ) = 2-F1 (a, b;c;w) 



r(c)r(A) ^ (a) w (6)„ _ 
r(a + A)r(6 + A)^ o n!(l-A)„ l J 



r( c ) 



r(a)r(6) 



to 



i) A E (a t^ w f^, A)w (i-^"[i»^-"-) 



n=0 



n!(n + A)! 



so that 



12. 



+ tp{a + A + n) + ip{b + A + n) - ^(n + 1) - ip(n + A + 1)] , (3.24) 

r( c )r(A) 



A(u,s) 



T(a + A)r(6 + A) 



^ s ) = r ( IlVX 1 [ ^ + A)+# + A)1 



(3.25) 
(3.26) 



11 Here and in related formulae, Tp(z) represents the Digamma function, denned by : ip{z) = yj^j- 
12 The terms tp(l) and ^>(A + 1) have been dropped from B because after dividing by A they amount to 
a contact term. 
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0„ 



("I) 



A+l 



r 2 (A) 

Putting this into Equation 3.17 we get: 



(o)a(6)a bl>(a + A) + ^(6 + A)] e 



A-l 



(3.27) 



G R (w,s) =4C S 



1 ( — 1) 



A+l 



r 2 (A) 



-(o)A(6)A[^(a + A)+V(6 + A)] . 



(3.28) 



3.3.2 Vector Correlators 

The calculation of the exact vector correlators in the K = — 1 case is similar to the scalar 
case: We write the expressions for the correlators in Equations 3.7 and 3.11 in terms of 
the variables defined in Appendix E and calculate the limits using the exact solutions. A 
detailed derivation is given in Appendix F, and here we quote the final results. 
For the longitudinal vector mode, we define: 



A = 



d-5 



(3.29) 



For the case of even d (A is non-integer), we obtain the expression: 



r R (l , o _ 4q rt 3 n2 r(-A)r(a + A)r(6 + A) 

a{,) ~d-5R d - 2qs T(A) T(a) T(b) 



(3.30) 



For the case of odd d (A is an integer), we obtain the expression: 



G^,s) = 2C v -^q : 



„d-3 
+ J2 



•1) 



A+l 



R d-2** r2(A + i 



-(o)a(6)a[^(o + A)+V(6 + A)] . 



(3.31) 



The other components of the longitudinal mode correlator can be calculated from Equa- 
tions 3.8 and 3.9. 

For the transverse vector mode, we define: 



A 



d-3 



(3.32) 



The caluclation gives for even d: 



r d-3 + A ) + A ) 



Gf ± (oj,v) = 4(d-S)C v -^ 



R d-2 F ( A ) r ( c 



r(6) 



(3.33) 



and for odd d: 



Gf ± (u,v) = 8C t 



d-3 



\A+1 



R d-2 F 2( A ) 



(o)a(6)a [ip(a + A) + ^(b + A)} 



(3.34) 
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3.4 Asymptotic Correlator Expressions (for K > — 1 ) 

Continuing Subsection 2.4, the goal of this subsection is to find an approximate expression 
for the correlators of the CFT operators coupled to the scalar and vector fields in the bulk 
spacetime for large values of the frequency, where |A| 3> q (corresponding to -Jr- = fixed 
and ji 3> 1) and ImA < 0. The general stages of the calculation are as follows: 

1. We start from the effective potentials and their poles calculated in Subsection 2.4. 

2. We solve the approximate QNM equations around each pole using Bessel functions. 

3. Using the asymptotic limit of |A| 3> 1, we replace the Bessel functions with their 
asymptotic forms. 

4. We match the solutions on each region of the complex plane, and apply the boundary 
condition at the horizon, to get an expression for the solution near the boundary up 
to a multiplicative constant. 

5. We use the approximate solution near the boundary and the expressions for the 
retarded correlators from Subsection 3.2 to obtain an asymptotic expression for the 
correlators. 

3.4.1 Calculation of Asymptotic Correlators for the Scalar Case 

Proceeding from the definitions in Subsection 2.4, the expression in Equation 3.5 can be 
written in terms of eft and z (after dropping the contact term and using the fact that 
z « z — 1 for z — > 1): 



o«( U ,.)--2cA ,( - 1) '** w(i) 



R d ~ z d-2 



(3.35) 



We continue in a method similar to the one applied in Subsection 2.4: find asymptotic 
solutions around the poles of the effective potential and "match" the solutions 13 . As before, 
we assume that |A| S> 1 and |A| S> q. We also assume that A is in the 3rd quadrant of the 
complex plane (Re A < 0, ImA < 0). The solutions around z — > zq (z — > oo) are given by: 

P+{z - zq) = V2vrA(5 - z )J (A(z - So)) (3.36) 
P (z - zq) = V2vrA(z - z )Y (X(z - z )) . (3.37) 

Replacing the Bessel functions with their asymptotic forms gives for direction (1): 



^4 e *(- A5 o-f ) _|_ j 4_ e «(-Azo-^) e «A5 _|_ i(Ai +f ) _ ^4_ e *( Az o-f ) 



e- lX ' z , (3.38) 



13 While in the general case, where Im(ASo) 7^ 0, the anti-Stokes lines emanating from 5 = and zq 
don't align, the "matching" of the solutions in the region (2) is still possible since there are no Stokes 
lines between these two anti-Stokes lines in this region, so that the asymptotic solution remains the same 
between them. 
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and for direction (2): 
6 ps \a+S-™o+^ - 3A.e l ^> + ^ 



e iXz + 



A e i(ASb+5) -^_ e i(A2o- 5) 



- «A5 



From the boundary condition at the horizon (z = 0, 2 — > — oo) we have: 

<j) ps e" iAi . 

There are now two possible cases: 
1. (i is even. In this case the solution around z — > is approximately: 

cf> ps B + P + (z) + B_P_(*) = B + y/2-K\zJ 3 ±{\z) + B^\/2ir\SJ_j ± (\S) 



(3.39) 
(3.40) 



(3.41) 



where ji = d — 1 (so that y is non-integer). Replacing the Bessel functions with 
their asymptotic form, we have: 



B + e~ il3 + + B. 



-iP- 



-iXz 



(3.42) 



where (3± = J (1 ± ji). We now match the solutions on lines (1) and (2). We have 
for line (1): 



For the section (2) we have: 



A. 



A i(-xz +^) _ 3A _ e i(-\Zo+V = B + e-*P+ + B. 



-%A A 



-0- 



^ +e i(A2o+|) _ j4 _ e i(A5o-f ) = B+e iP+ + B _e ifi - . 
Putting Equation 3.43 into Equations 3.44 and 3.45 we get: 

B + e~ ie+ + B^e" w - = 2iA + 
B+e i9 + + B_e w - = 4iA + , 

where 9± = Xzq — f — j3± . Solving for B± we can get an expression for : 



(3.43) 

(3.44) 
(3.45) 

(3.46) 
(3.47) 



2 e- i9 ~ 
4 e ie ~ 



e- ie + 2 
e' 



J0+ 4 



J2i6- 



3 2i6L 



+ 2 



(3.48) 



(here we used the fact that d is even). 

Next we evaluate the correlator. Developing <p around z = we have: 
cj) PS B + V2tt\zJ 3jl (\z) + B^2n\zJ^r L (\z) 

2 2 

psP v^A.i'A:) ; ^ l "f A - ' 



(-lA 2 z 2 ) fe 



+ J B+v / 2vrAz(-Az)^ 



(3.49) 
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Defining the coefficients: 



A(u,s) 
B(u,s) 



2ttA- a+ I 



2- A r(-A + 1 



2 A r(A + 1) 



(3.50) 
(3.51) 



(where A = ^ = we S e * : 



(2A)^4^- 1 



(d-1) 



iXV' 1 T(-A) e : 



2i0_ 



d-2 



T(A) e 2ie - + 2 



(3.52) 



And finally putting this into Equation 3.35: 



rf 1 r(-A) A'A 



G K (w,s) « -2(d-l)C s 



R d T(A) V 2 



d-1 „2i0_ 



3 2i6L 



+ 2 



(3.53) 



(Note that this expression is true up to the addition of contact terms). This asymp- 
totic formula is only true for Re A < 0. In order to get an expression for Re A > 0, 
we can make use of the symmetry properties of the correlator 14 : 



G r (oj,s) = G r *(-uj*,s) 



We get for Re A > 0: 



G R (u,s) = -2(d-l)G 



R d T(A) V 2 



-1 e 2i 



+ 2 



where we define: 



7T 



.( A ) = -^(-A*) = A^- T (d-3) . 



(3.54) 



(3.55) 



(3.56) 



Looking at the asymptotic behaviour of the expressions for Re A < and Re A > 0, 
we may unify them into one asymptotic expression: 



G R (u,s) ~ 2{d-l)C t 



r d-i r(_ A ) fi\ 



-1 e 2i 



2 e 



2i0- 



R d T(A) \2 J e 2id - + 2 e 2i0_ + 2 



(3.57) 



Finally, using A 



wR 2 

r + 



, the asymptotic expression may be written: 



G r (uj,s) ps 2{d-l)C s R 



d _ 2 F(-A) / iu 



- 



r(A) V 2 



d-1 „2i0_ 



2 e : 



2i6- 



+ 2 e 2*0- + 2 



(3.58) 



2. d is odd. In this case the solution around z — > is approximately: 



B + P + (z) + £_P_(z) = B+v^vrAzJ^A-z) + £_ V^TrAlYn (AS) 

2 2 



(3.59) 



This property can also be deduced from the QNM equations themselves. 
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where j\ = d — 1 (so that y is an integer) . Replacing the Bessel functions with their 
asymptotic form, we have: 

4> « \(B + - iS_) e - i/3 + 



e i\z + 



(B+ + iB^)e^+ 



-i\z 



(3.60) 



where /3+ = ^(1 = f<2. We again match the solutions on lines (1) and (2). For 
(1) we again have: 



For (2): 



A i(_Azo-f) + A _j(-^o-^) = o _ 



A+e i(-Xzo+'f) _ 3A _ e i(-\z + l ) = (_g + _ 

^ +e *(Aio+f ) _ j4 _ e i(A5 -f ) = ( 5+ + iB _y 



(3.61) 

(3.62) 
(3.63) 



Putting Equation 3.61 into Equations 3.62 and 3.63 we get (defining 9 + = \z$ — j 



i.e. 



(B+ - iB-)e 
(B+ + iB)e- ie + = 2iA+ . 



(3.64) 
(3.65) 



Solving for B± we can get an expression for 



4 -ie ie ^ 
2 ie~ ie + 



e w + 4 
e~ ie + 2 



2z- 



(3.66) 



Next we evaluate the correlator. Developing <j> around z = 0we have (we again define 
A = 4f = ^r 1 ): 



S + \/2^AIJa(Az) + 5_v / 2^AlY A (Ai 



Defining the coefficients: 



^(w,s) 



vr2- 



-A+l 

3— (A-1)!S_ 



(3.67) 



(3.68) 



^' S) = vr2Ar(A + l) ln 



we again have: 



2^A" 



2^A A+ 2 



^(A)! ^ 1 ) + ^ + 1 )]^- + 2 A r( A + i) 



B A 



d- z (, 



(2A)*4e 2A - 1 
v J A(u>,s) 



(3.69) 



(3.70) 
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And finally after putting this into Equation 3.35, and again making use of the sym- 
metry properties of the correlator we obtain: 



G r (uj, s) « 4C S R 



1-2 ' 



-1 



,A+1 



T 2 (A) V 2 



d-1 



,2i0 4 



2vu- 



,2i0+ 



e 2ie + , iX 

2m — = h 2 In — 

e 2ie + _ 2 1 2 



(3.71) 



where: 



7T 



(9+(A) = -6>;(-A*) = Xz* + -(d+ 1) . (3.72) 
3.4.2 Calculation of Asymptotic Correlators for the Vector Case 



The calculation of asymptotic expressions for the vector correlators is similar to the scalar 
case. A detailed derivation is given in Appendix H, and here we quote the final results. 
For the longitudinal vector mode, we define: 

d ~*. (3.73) 



A = 

2 

We distinguish between two possible cases 15 : 

1. d is even. In this case we obtain the asymptotic expression: 



GS(",s) 



4C v ^-4/ r FI\2 r (- A ) f iuJ 



d-5 



R«-\L{ 



r(A) 

c 



d-5 



2,H - ■ 2 cos I f., I 



2i6- 



+ 2 COS 



d-2 



^2i8- 



2 COS | I < 



2i6- 



2 cos 



d-2 



where: 



7T 



= Xz + -(d-7) 



TT 



e„(x) = -et(-x*) = Xz* --(d-7) . 

2. d is odd. In this case we obtain the asymptotic expression: 



(3.74) 

(3.75) 
(3.76) 



d-5 



2vri- 



3 2i0 4 



+ 2 cos 



d-2 



2m- 



-,2id i 



+ 2 cos 



d-2 



+ 2 In 



iA 



(3.77) 



lj Note that for the case of vector gauge field perturbations with d — 4, the effective potential for both 
modes is V(z) = q 2 g(z), and has a pole of degree — | at 5 — > zq. In this case the method applied here 
can't be used since the equation is not of the Bessel type around the pole So, and since q obviously can't be 
neglected (neglecting it results in an an expression with no poles). This case therefore requires a separate 
and more complete treatment, which shall not be investigated here. This section will therefore focus on the 
case of d > 5. 
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where: 



7T 



?+ = Xzq — — (d — 3) 



7T 



?+ (A) = -^(-A*) = A^ + i (d-3) . 



(3.78) 
(3.79) 



The other components of the longitudinal mode correlator can be calculated from Equa- 
tions 3.8 and 3.9. 

For the transverse vector mode, we define: 

d — 3 , 
— — . (3.80) 



A 

2 

We again distinguish between two possible cases: 

1. d is even. In this case we obtain the asymptotic expression: 



Gxi(w,v) 



4(d - 3)C V R 



d _ 4 T(-A) / iuj 



r(A) 



d-3 



2 COS | I 



2i0_ 



2 cos 



d-2 



e W- + 2 cos ( ^ e 2ie - + 2 cos , (/ _., 



where: 



Xz 



7T 



(d-5) 



?_(A) 



7T 



r ( _A*) = A^- T (d-5) . 



(3.81) 

(3.82) 
(3.83) 



2. d is odd. In this case we obtain the asymptotic expression: 




where: 



= Az 



7T 



(d-1) 



7T 



? + (A) = -«;(-A*) = A^ + -(d-l) . 



(3.84) 

(3.85) 
(3.86) 



Notice that, taking into account the differences in the definition of A and 9± between 
the longitudinal and the transverse perturbation modes: 

. long 



A 



^ long 



a trans 
^trans 



- 1 
7T 



(3.87) 
(3.88) 



the asymptotic expressions for G^j_(u;,v) and G^(oj,s) are identical. This is expected 
and serves as a check for our calculations, since in this limit we have neglected the spatial 
mode L^T, so that the differences between the longitudinal and the transverse modes are 
negligible. 
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3.5 Numerical Calculation of Correlators 

In the following we present some results of exact numerical calculations of the gauge theory 
retarded correlation functions for several cases, and compare them to the asymptotic an- 
alytical expressions (as discussed in Subsection 3.4). The correlation functions have been 
calculated using the method outlined in Appendix I. 

The following calculations were made for a degree of N = 200 and zq = 0.6 (see 
Appendix I for details): 

• The case of scalar perturbation in d = 4 bulk dimensions, with K = 1 (spherical 
boundary topology) and q s = (s-wave). The results are given in Figures 8a, 8b, 
8c and 8d. 

• The case of longitudinal vector perturbations in d = 5 bulk dimensions, with K = 
(flat boundary topology) and q s = 0.6. The results are given in Figures 9a, 9b, 9c 
and 9d. 

• The case of transverse vector perturbations in d = 5 bulk dimensions, with K = —0.5 
(hyperbolic boundary topology) and q v = 1. The results are given in Figures 10a, 
10b, 10c and lOd. 

For the sake of better visual representation of the results, we define the "normalized" 
retarded correlators 1(A) by "normalizing out" the power and logarithmic behaviour of 
the correlators, keeping only the part of the correlators that behaves as O(l) for |A| —> 
oo 16 . These "normalized" correlators have poles at the QNM frequencies, and approach a 
constant on each side of the QNM frequencies line for frequencies far enough away from 
that line. 

Figures 8a, 9a and 10a show the value of |/(Ao + sAA)| as a function of s, where Ao 
and AA are the analytical asymptotic parameters of the QNM spectrum as calculated in 
Subsection 2.4 and s is a positive, real parameter. As expected, we find the poles of the 
correlation function approximately at integer values of s. 

Figures 8b, 9b and 10b show the value of Re / (A^e - * 9 ) as a function of 9, where A# 
is some constant radius in the frequency complex plane (in this case Xr = 20). As ex- 
pected, the numerically calculated values are close to the ones of the analytical asymptotic 
expressions calculated in Subsection 3.4. 

Figures 8c, 9c and 10c show a contour plot of Re 1(A) in the 4-th quadrant of the 
complex frequency plane. Both the poles at the QNM frequencies and the step-function- 
like behaviour of the "normalized" correlation function can be clearly seen in the plots. 

Figures 8d, 9d and lOd show the values of the appropriate spectral functions on the 
real frequency axis (p(X) in the case of scalar perturbations, ptt(X) in the case of longitu- 
dinal vector perturbations and p±±(X) in the case of transverse vector perturbations). At 
large frequency values, these functions approach the analytical asymptotic expressions and 
are proportional to X s (where s is the appropriate exponent for the perturbation type), as 

16 See Appendix I for a more precise definition. 
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IKAo+sAA)| 



Scalar d=4 K=1.00 q=0. 



Scalar d=4 K=1.00 q=0. |A|=20.0 



ReI(]A|Exp[-i9D 



23 43 



63 83 




- en 



(a) 



(b) 



Scalar d=4 K=1.00 q=0. 




Scalar d=4 K=1.00 q=0. 




(d) 

Figure 8. Numerically calculated "normalized" retarded correlation function for the case of scalar 
perturbations with d = 4,K — 1 and q s = 0. The plots show: (a) The absolute value of the "nor- 
malized" correlation function along the analytically calculated asymptotic line of QNM frequencies, 
given by the parameters: A = -2.7206 - 2.7205i AA = -1.9691 - 2.3502vi, (b) the real parts of 
both the numerically calculated and the analytical asymptotic "normalized" correlation functions 
on a circle in the complex plane of radius |A| = 20, (c) a contour plot of the real part of the nu- 
merically calculated "normalized" correlation function in the 4-th quadrant of the complex plane, 

r d ~ 1 

(d) the spectral function associated with G R along the real frequency axis, in units of — C s -^a-. 



expected from conformal symmetry considerations. In the longitudinal vector case (Fig- 
ure 9d), the spectral function oscillates at lower frequencies, due to the existence of the 
hydrodynamic diffusion pole. 
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(b) 



LongVector d=5 K=0. q=0.600 




LongVector d=5 K=0. q=0.600 




(d) 

Figure 9. Numerically calculated "normalized" retarded correlation function for the case of lon- 
gitudinal vector perturbations with d = 5,K = and q s = 0.6. The plots show: (a) The absolute 
value of the "normalized" correlation function along the analytically calculated asymptotic line of 
QNM frequencies, given by the parameters: Ao = AA = — 2 — 2i, (b) the real parts of both 
the numerically calculated and the analytical asymptotic "normalized" correlation functions on a 
circle in the complex plane of radius |A| = 20, (c) a contour plot of the real part of the numerically 
calculated "normalized" correlation function in the 4-th quadrant of the complex plane, (d) the 

spectral function associated with Gf t along the real frequency axis, in units of — C, 



V R d- 
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(b) 



TransVector d=5 K=-0.500 q=1.00 




TransVector d=5 K=-0.500 q=1.00 




(d) 

Figure 10. Numerically calculated "normalized" retarded correlation function for the case of 
transverse vector perturbations with d — 5,K — —0.5 and q v = 1. The plots show: (a) The absolute 
value of the "normalized" correlation function along the analytically calculated asymptotic line of 
QNM frequencies, given by the parameters: A = -1.4142 - 2.0CXXK AA = -1.4142 - 2.0000i, 
(b) the real parts of both the numerically calculated and the analytical asymptotic "normalized" 
correlation functions on a circle in the complex plane of radius |A| = 20, (c) a contour plot of the 
real part of the numerically calculated "normalized" correlation function in the 4-th quadrant of 
the complex plane, (d) the spectral function associated with G^ ± along the real frequency axis, in 

r d-3 

units of -C„ jj3=5 ■ 
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4 Discussion 



In this work we studied some of the properties of the retarded correlation function of 
composite operators of finite temperature gauge field theories with dual gravity descriptions 
in the large coupling limit. Of specific interest was the dependence of these properties on 
the topology of the space on which the gauge theory is defined, and its temperature. This 
goal was achieved by performing exact analytical calculations for specific cases where it was 
possible (the hyperbolic case with T = T c ), by performing numerical calculations and by 
obtaining approximate expressions for the large frequency limit {—jr- = fixed and j>1). 

In the large frequency limit we saw that in all cases (all topologies, temperatures and 
perturbation types) the correlation functions exhibit the same kind of asymptotic pole 
structure, with a constant "gap" between two consecutive poles. Moreover, we saw that 
asymptotically these poles are produced by Bose-Einstein-like factors in the correlation 
functions, with some complex parameter multiplying the frequency. 

We have also seen that the spatial topology on which the gauge theory is defined is 
indeed noticeable in the properties of the correlation function and its poles (corresponding 

to quasinormal modes in the dual bulk spacetime), in the case where v jp — = fixed. In 
particular, the asymptotic "slope" of the correlator poles encodes information about the 
spatial topology and the temperature, but is "universal" otherwise - it doesn't depend 
on the type of operators considered. In the spherical (positive curvature) and flat (zero 
curvature) cases, this slope always has a real part and an imaginary part. In the spherical 
case, it depends on the relation between the spatial curvature and the temperature, while 
in the flat case it's constant (as expected, since in this case the temperature is the only scale 
in the theory, and the correlator poles are necessarily proportional to it). In the hyperbolic 
(negative curvature) case, the slope has a real part for temperatures above T c , which goes 
to in a non-smooth way for T — > T c (In d = 4, 5 dimensions the real part is also for 
temperatures below T c ). At T = T c the slope is completely imaginary, and exact analytical 
expressions can be found for the correlation functions and their poles. Interestingly, the 
poles "gap" in this case, and indeed the general form of the correlation functions, is exactly 
what one would expect to find for the retarded correlation function of composite operators 
in a free thermal field theory (the poles "gap" in this case corresponds to the Matsubara 
frequency "gap" - see [24]). 

Several open issues are left for further investigation: 

• This work has dealt with uncharged bulk black hole solutions, but we can also expand 
the scope and discuss charged black holes, that correspond to gauge theories with 
non-zero chemical potentials (see [15]). The background of the perturbations in such 
solutions contains both a black hole metric and a non-zero gauge field. We can then 
ask - how does the charge change the asymptotic properties of the gauge theory 
correlation functions and their dependence on the spatial topology? 

• As noted in Subsubsection 3.4.2, the method we have employed for calculating the 
asymptotic expressions for the correlation functions doesn't work as is for the vector 
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perturbation modes in d = 4 bulk dimensions. This case requires a different treat- 
ment using complex WKB analysis, although numerical calculations show that our 
conclusions regarding the asymptotic "slope" of the correlator poles remain true in 
this case as well. 

• While the asymptotic structure of the correlators (and their poles) has a geometrical 
meaning in the gravity (bulk) side, it is unclear if it has an interpretation in the 
strongly coupled dual field theory. Weak coupling calculations don't seem to exhibit 
such a structure (see [24]). The asymptotic poles "gap" in the large coupling limit 
depends only on the spatial topology and the temperature, suggesting it may have an 
interpretation involving them, but it can't be trivially deduced from the symmetries 
of the theory alone. 

• Considering the properties of the correlation functions in the hyperbolic case for 
T — > T c , does this temperature have any importance in the dual gauge theory? 
In [26], it is suggested that in d = 4, the behaviour of the QNM frequencies slope 
near this temperature is related to the known phase transition between the topological 
(hyperbolic) black hole and the MTZ solution (a hyperbolic black hole with scalar 
hair - see [25], [27]). It is not clear, however, what is the interpretation of such 
a phase transition in the gauge theory side, and whether the argument extends to 
higher dimensions. 
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A List of Notations and Definitions 

Here we list for convenience some of the main notations used throughout this paper. 

• d is the number of bulk dimensions. 

• R is the AdS radius. 

• Rq = R 2 R^ T is the scalar curvature of the Qd-2 manifold. 

• r + is the radius of the bulk black hole horizon. 



• oj is the frequency of a perturbation mode, L^T = is the Laplace operator 

eigenvalue of the perturbation mode on the dual field theory side, for scalar or vector 
respectively. 
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• t c - = h\jjMU ( for the case of k < °)- 

• = 1 + K(1- zf -(1 + K)(1- zf- 1 . 

• z = - J2k=l 7fe ln (! - ^fe)> where g{z k ) = and ^ k 



g'(zk) 



for A; = 1, . . . , d — 1. 



C s /„ is the normalization constant for the bulk action of the scalar/vector field re- 
spectively. 

scalar perturbation, 
transverse vector perturbation, 
longitudinal vector perturbation. 




"(d-5)-f ±i^2_i (d _ 5) 

^s 2 -i(d-3) 



2 31 2 

iX _i_ 
2 31 2' 



scalar perturbation, 
transverse vector perturbation, 
longitudinal vector perturbation. 



A^o + f — 3) scalar perturbation, 

A^o + f (d — 5) transverse vector perturbation, 

A^o + \ (d — 7) longitudinal vector perturbation. 

Aig — | (d — 3) scalar perturbation, 
\2q — \ (d — 5) transverse vector perturbation, 
A5q — f (d — 7) longitudinal vector perturbation. 



0-(A)H 



MA) = -r (-a*) = < 



Azo — f (d + 1) scalar perturbation, 
0+(A) = ^ Azo — | (d — 1) transverse vector perturbation, 
Xzo — \{d — 3) longitudinal vector perturbation. 

Azq + f (d + 1) scalar perturbation, 



"(A) = -« 



-A*: 



A5q + j(d — 1) transverse vector perturbation, 
Xzq + j(d — 3) longitudinal vector perturbation. 



B Thermodynamic Quantities 

Given the black hole metric, one may calculate the usual thermodynamic quantities re- 
lated to the black hole. In the context of AdS/CFT, these quantities are related to the 
respective quantities in the field theory side. Defining: p = '-^ and K = we have for 
the temperature: 

f(r + ) 



4tt 



4ttR 



[2 + (K + l)(d-3)} 



(B.l) 



For the total energy we have: 



E = EFT JAz^m pd - 1{K+1)+Eo , 



(B.2) 



WirGR 

where Vol(fl^ 2 ) 1S the total volume of the ^5T 2 manifold, G is the bulk gravitational 
constant and Eq is independent of the black hole parameters. Finally the total entropy is 
given by: 



VoWd-2) d -2 d-2 
D = JFT = r + = TF< P 



(B.3) 



AG + AG 

For the sake of discussing the field theory correlators and the quasinormal modes, we shall 
also define the following dimensionless quantities, for an oscillation mode of frequency uj 
and Laplace operator eigenvalue L s : 

ujR d — 1 uj I 

P 



A 



2tt T 



1 + 1/1 



d-l 



L, L^ T R d-lL 



FT 



P 



2vr T 



1 + 1 1 



d-2 (2nT) 2 
1 

d-l K n 
d-2 (2ttT) 2 



(B.4) 
(B.5) 



Table 1 contains a summary of the physical quantities defined in the bulk spacetime, the 
corresponding field theory quantities and their dimensions. 

TR(p) for d=5 





Figure 11. TR as a function of p for the spherical (k = 1), flat (k — 0) and hyperbolic(fc = — 1) 
cases in d = 5 bulk dimensions. 

C QNM Equations Derivation 
C.l QNM Scalar Equation 

The EOM for a minimally coupled, massless scalar field $ are given by 

V M V$ = , (C.l) 
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Table 1. A summary of physical quantities in the bulk and field theory spacetimes and their 
dimensions 



Bulk Quantity 


Field Theory Quantity 


Quantity 


Dimension 


Quantity 


Dimension 


t 


[L] 


t 


[L] 


x l 


[L]° 


oo — Rtjc 


[L] 


r 


[L] 






ds 2 


[L? 


dSprp 


[L? 


R 


[L] 






r+ 


[L] 






p = r+R- 1 


[L}° 






Rn 


[L]° 


Rq T = RqR 2 




k 


[L}° 


k FT = kR- 2 


[L]- 2 


K = kp- 2 


[L}° 






T 


[L]- 1 


T 


[L]- 1 




[L]- 1 






L s 


[L}° 


Lf T = LgR^ 1 


[L]- 1 


A = uRp- 1 


[L}° 






q s = Lsp' 1 


[L}° 






G 


[L] d - 2 






E 


[L]- 1 


E 


[L]- 1 


S 


[L}° 


S 


[L]° 



or 



-fJ 



A (V^gg^d^) = o 



(C.2) 



(where g = det(g^) ). 

To find a QNM, we look for solutions of the form 

§(M,x)=^(z)e-%(x), (C.3) 

where H L 2(x) is a (scalar) eigenfunction of the Laplace operator defined on the manifold 
with corresponding eigenvalue of —L 2 . Putting this in, and noticing that 



(C4) 



(where d and 5 are the exterior derivative and codifferential defined on 0^_2, respectively), 
we obtain Equation 2.9 (scalar QNM equation). 
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The QNMs in asymptotically AdS spacetimes are denned as the solutions to the equa- 
tions of motion obeying the incoming wave boundary condition at the horizon and the 
Dirichlet boundary condition at the AdS boundary (at least in the context of AdS/CFT - 
see Subsection 2.1 and also [13], [17], [22]). In terms of the defined variables and coordi- 
nates, these translate to 

i>\ z=0 ~z~% VU = 0, (C.5) 

as in Equation 2.10. 

C.2 QNM Vector Equations 

The EOM for a vector (gauge) field A are given by 

V M F^ = , (C.6) 

or 

-q^ (y=^V^/?) = . (C.7) 

Putting the metric from Equation 2.5 Equation C.7 we get for the (3 component of the 
equation: 

+ i E -7=<* (v«/ 7 ^7) = • (C8) 



r t? Van 

To find a QNM, we look for solutions proportional to eigenfunctions of the Laplace 
operator. The t and z components would be proportional to the eigenfunctions of the scalar 
Laplace operator: At,A z oc H L 2(x), where H L 2 is defined as in Equation C.4. The other 
components would be proportional to the corresponding components of the eigenfunction 
of the vector Laplace operator: A, oc A £ 2 j(x), where A L 2 is defined as follows: 

A n A i2 (x) = 5dA L2 (x) = ]T -^=d 3 (VanF%) = L^ i2 (x) . (C.9) 



There are now two possibilities that correspond to two different modes: 

1. The "longitudinal" mode. This is the case where L v = , or 5dAo = 0, from which 
we can deduce 

( dA , dA ) = (A , 5 dA ) = => F = dA = . (CIO) 

Ao is, therefore, a closed form and can be written as Ao,i = di<f), where <f> is some 
(scalar) function 17 . To get the equation for a QNM we choose <j> to be an eigenfunction 



17 More precisely, Aq can be written as a sum of an exact form d(f> and some harmonic form B. The 
harmonic component satisfies <5B = 0, and can therefore be treated as a "transverse" mode. 
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of the scalar Laplace operator <j> = H L 2 . The solution will then be of the form 

/ A t (t,z,x) \ ( f(t,z)H L 2(x) \ 

h(t,z)H L2 (x) 



A z (t,z,x) 
Ai(t,z,x) 



k(t,z)diH L2 (x) 



(C.ll) 



\A d _ 2 (t,z,x.)J \k(t,z)d d - 2 H L 2{x)J 

Putting this into Equation C.8 we get the following equations: 
For the t component- 

1 

O^z) 



I = g(z)(l - z) d - 4 d; 
For the z component- 

U z 



W-4 ( 9 * h ~ 9 *f) 



q 2 s (d t k - /) = . (C.12) 



ri 1 



p 4 g{z) 



d t {d t h - d-j) - g 2 s (h - d z k) = o 



(C.13) 



From Equations C.12 and C.13 we can take d z (I) + dt(II) = 0, and defining 

(dth — d z f) = -0(z)e _iw ', we finally obtain Equation 2.11 (longitudinal vector QNM 

equation) . 

As for the boundary conditions, the incoming wave boundary condition at the horizon 
implies that 

VL=o ~ , (C14) 
while the Dirichlet boundary condition at the boundary implies 



f,h,k\ z=1 = 



(C.15) 



Setting z = 1 in Equation C.12 and applying the boundary condition we get for ip 
the condition: 



(i - zf-^d. 



i 



(l-z) 



d-4 



. 



(C.16) 



2=1 



2. The "transverse" mode. In this case L v ^ 0. Applying the codifferential 5 to Equa- 
tion C.9 we get 

(C.17) 
(C.18) 



0=55 dA L 2 = L 2 V 5A L 2 , 

and since L v / we see that 

5A L2 = -J2-^di{V9nAi 2 ) 



fJQ 

The solution in this case will be of the form 
/ A t (t,z,x) \ ( 



. 



A z (t,z,x) 








a(t, z)A L 2 1 (x) 



(C.19) 



\A d _ 2 (t,z,x.)J \a(t,z)A L 2 td _ 2 (x)j 
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Putting this into Equation C.8, the equations for the t and z components are auto- 
matically satisfied accroding to Equation C.18. 
Taking into account Equation C.9, we get for the i component 



(1 - z) d ~ 4 d, 



9{z) - \ ri 1 
d z a 



(1-z) 



p 4 g(z. 



d?a - = . (C.20) 



Setting a(t,z) = e lU)t ip(z) we finally obtain Equation 2.15 (transverse vector QNM 
equation) . 

The boundary conditions are again given by the incoming wave boundary condition 
at the horizon and the Dirichlet boundary condition at the boundary, so that: 



(C.21) 



D Derivation of General Formulae For Correlators 
D.l Scalar Correlators 

Consider a massless, minimally coupled scalar field in the bulk spacetime (where the bulk 
metric is given by Equation 2.5). The action of the scalar field is given by: 



S[<f>] = C S dz dt dx 1 ... dx d - 2 ^gV^ 

JO J — OO 



(D.l) 



where C s is the appropriate normalization constant for the bulk scalar field. Since (ft is 
on-shell, it satisfies the EOM: 

v^vv = o , 



(D.2) 



so that, using integration by parts: 
S = C S 



pi poo p 

dz dt dx 1 ... dx d ~ 2 ^g (V M {(ftV^cft) - (ftV^V^ 

JO J-OO 

poo p 

C, I dt dx 1 ... dx d ~ 2 ^~gg z 

-DO J^d~2 



Z = l 



2=0 



(D.3) 



Using the metric from Equation 2.5 gives: 



„d-l 



so that 



(l-z) d 

dt I dtt d _ 2 (, 



z=l 



z=0 



(D.4) 



(D.5) 



(where d^d-2 = y/gadx 1 ■ ■ ■ dx ). 

Let H s {x) be a complete basis of eigenfunctions of the Laplace operator defined on the 
£ld-2 manifold, where s here stands for all of the indices (discrete or continuous) required 
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to uniquely identify the eigenfunctions. The normalization of these functions will be chosen 
so that: 

/ dn d _ 2 H* s ,(x)H s (x) = 6 SS , , (D.6) 

where 5 ss i here stands for the Kronecker delta in case of a discrete spectrum and a Dirac 
delta in case of a continuous spectrum. One may then expand the field (ft in terms of these 
functions and the Fourier components so that: 

(ft(z, t,x) = f ^ e- luJt H s (x)^, s (z) , (D.7) 

J s 

where Yl s stands here for either a sum in case of a discrete spectrum or an integral (with 
the appropriate measure) in case of a continuous spectrum. Since (ft is a solution for the 
EOM, as stated in Subsection 3.1, ip u>a would be a solution of the scalar QNM Equation 
(Equation 2.9), with the incoming- wave boundary condition at the horizon: 

V^,s| 2=0 ~ z~ % c . (D.8) 

As for the boundary, notice that the metric of CFT spacetime is actually given by Equa- 
tion 3.4, so that in the CFT spacetime the normalized eigenfunctions H^ T {x) are given 
by 



-2 



I dn^ 2 H^*(x)Hr(x) = R d - 2 j <M d „ 2 H™{x)H! T {x) = 5 a 

=¥ H^ T (x) = R~^H s (x) . (D.9) 
As a consequence, on the boundary the field (ft may be written as: 

<ft(l,t,x) = f ^^e-^-^if s (x)^ ^,s) (D.10) 

J s 

(where (fto(oj, s) are the components of (ft in the CFT spacetime), which means the boundary 
condition on ip u>a at the AdS boundary is: 

^w,s| z _>i = iT^0 o (u;,s) . (D.ll) 

Defining a normalized i}) wa by tpus = R ^~ < fto( U} : sjV'w s> the normalized function satisfies 
the boundary conditions: 



~ Z c lftu,s 

z=0 



Putting Equation D.7 into the expression for the action, we have: 



= 1 . (D.12) 

z— >1 



3 i ^'-^ t H* s ,(x)H s (x)^ s '(z)d z ^ s (z) Z 1 . (D.13) 

z=0 
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Using 

J dte^'-^ = 2k8{J -oj) 
and the orthogonality of H s (x), the action becomes: 

2=1 



(D.14) 



2=0 



r d ~ l g(z)(l-z) 2 1 [ dw 
(1 - z) d «2 ^2 y 



E 



0o( w ) s )C, s (^)^^,s(^)</'o(w,s) , (D.15) 



z=l 



2=0 



so that 



where 



5[0] 



2=1 



2 = 



r^ 1 0(2) - - 



R d (1-z) 



Using the prescription described in [13], the retarded correlator is then given by: 

„d-l 



G R (to,s) = -%F(u,a,z)\ s ^ x = -2C S 



D.2 Vector Correlators 



(D.16) 



(D.17) 



(D.18) 



Next consider a vector gauge field in the bulk spacetime. The action of the gauge field is 
given by: 

fl poo f 

(D.19) 



/*! fOO r 

S[A] = C V dz dt dx 1 ... dx^^gF^F^ , 

JO J— OO 

where C v is the appropriate normalization constant for the bulk vector field. Since is 
on-shell, it satisfies the EOM: 

(D.20) 
(D.21) 



V^F^ = 



so that 



F„ V F^ = 2V M (A„F^) 



and using integration by parts: 



S — 



dt 



dx 1 ... dx d ~ 2 ^~gA v F z 



2=1 



2=0 



(D.22) 



Using the metric from Equation 2.5 gives: 



S = 2C V + g zz I dt I dn d „ 2 [g a A t F zt + g^ AF zj ] 



z=l 



z=0 



(D.23) 
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A z (z,t,x) = [ ^^ e -^ s (x)V s (^) (D.28) 



(where dfi rf _ 2 = ^/gn dx 1 . . . dx d 2 ). 

Let H s (x) again be a complete basis of eigenfunctions of the Laplace operator on the 
£ld-2 manifold, normalized so that: 

J dn d „ 2 H*,(x)H s (x) = 5 SS , . (D.24) 

Let us further define A v ( X ) clS cl set of vector eigenfunctions of the Laplace operator, such 
that: 

SdA v = LlA v L v /0, (D.25) 

and 5A V = 0. Here v stands for all of the indices required to uniquely identify the 
eigenfunctions. The normalization of these functions will be chosen so that: 

J 6n d - 2 g^ ti (x)A^ d (x) = <W . (D.26) 

~ 0-ff (x) 

Note that A„ along with - ^ ' form a complete orthonormal basis. One may expand the 
field A in terms of these functions and the Fourier components: 

A t (z, t,x) = I ^ J2 e-^HsWf^iz) (D.27) 

du 
2^T 

s 

^ s J v 

Since A is a solution for the EOM, as stated in Subsection 3.1, f w<s , /i WjS) &w,s and a WjV 
would be solutions of the vector QNM Equations (see Appendix C) with incoming-wave 
boundary conditions at the horizon. As for the boundary, since the field theory metric is 
given by Equation 3.4, in the field theory side the normalized eigenfunctions are given by: 

H^ T (x) = R-^H s {x) (D.30) 

A^ T (x) = R~^A V . (D.31) 

As a consequence, on the boundary the field A may be written as: 

A t (l,t,x) = I ^^e-^R-^HsCaOA^s) (D.32) 

^ s 

*<!■«.*)= 

J s s 

+ / ^£e-^-^i v ,(*K(u,,v) (D.33) 

(where A®(uj, s), -A{j(w, s) and -Aj_(u;, v) are the components of A in the field theory space- 
time). The boundary conditions at the AdS boundary are therefore: 

f^s\ z ^ 1 = R- d -^A° t (oj,s) (D.34) 

LL, = R-^L^AVAu, s) (D.35) 



Ow.vL-n = R-^A° ± (uj,v) . (D.36) 
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Next, one may derive expansions for the Fn V tensor: 

F zt = - [ ^ Yl e ""^(*) [-iuK* ~ d t f Ufi ] (D.37) 

J s 

F Z j = [ e-^djHsix) [d z k^ s - h U;S ] + f^J2 e- iuJt A vJ (x)d zauliV . (D.38) 

J s ^ v 

We proceed by inserting these expressions into the expression for the action. Using the 
orthogonality properties of Equations D.24 and D.26, and also the following derived prop- 
erties: 



tn d -29%diH${x)djH a {x) = (dH s ,, dH s ) = (H s ,, 5 dH s ) = L Z J S 



(D.39) 



dO d _ 2 g%aiH}(x)A V j(x) = ( dH B ,,Av) = (H s ,, 5A v j = . (D.40) 

We see that the action decomposes into a longitudinal part S», that involves / WjS , h UjS and 
kuj^, and a transverse part S± that involves only a w v , so that: 



S = Sn + S± . 

For the longitudinal part of the action, we get the following expression: 



(D.41) 



(i - z y 



[1 - z? t * 



p 2 g{z) 



+ 



R 2 



9s ^lu s (&zkiu,s ^u;,s) 



z=l 



z=0 



(D.42) 



(where the definition q s = was used). Since / w>s , h w>a and k^ tS satisfy the longitudinal 
vector QNM Equations, we may use Equation C.13 to get the relation: 

^2 



Using this relation, along with the definitions: 

~X,U1,S = 2<-<-'^CJ,S fu},S ) 

we obtain the expression: 



„d-3 



S\\ — 2 Cy 



(1-2) 



d-4 



z=l 



z=0 



(D.43) 



(D.44) 
(D.45) 



(D.46) 



Looking at Equations C.12 and 2.11, we see that Xco,s an d tpw,s satisfy the following system 
of differential equations: 



~g(z)(l-z) d - A d, 



1 



(1 - z) d " 4 



,2- 



A 2 - qig(z) 



dzXu 



(D.47) 
(D.48) 
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From Equations D.35 and D.34 we can derive the boundary condition for Xw,s : 

Xu,s\ z ^i = -wR'^L- 1 Af ] (u>,s)-R-^A t (uj,s) . 
Next, define the normalized x^,s and ij) u s by: 



-iuR- — L- 1 A? l (u},s) - R- — A° t (io,s) 



iuR-^L^A^u, s) - R-^AKu, s) 



Xoj,s 



(D.49) 

(D.50) 
(D.51) 



so that Xlo,s an d ip u,s also satisfy Equations D.47 and D.48, with the boundary condition: 

X Wl sUi = 1 • (D.52) 
Using these definitions in the expression for the action gives: 
do; 



Si 



2tt 



A 



i—A?,(u),B) + A$(u,B) 



A 



where: 



„d-3 



J"||(a),s, 2r) = — 2C„- 



2 = 1 



z=0 



(D.53) 



(D.54) 



R d - 2 (l-z) 

Using the prescription described in [13], we can then get the components of the retarded 
correlators that correspond to the longitudinal mode: 



Gtt(w,s) = -2Tn(ui,s,z] 



r d-3 -, 



„<2-3 



4C„ 



1 



R d-2 _ ^ 



(D.55) 



G*( W ,s)=i-G£(u,,s) 



Note that from Equation D.47 we have: 



(1 - z) d ~ A d z 



(l-z) 



.1 -4 YUJ,s 



Z-+1 



(D.56) 
(D.57) 



(D.58) 



Turning to the transverse part of the action, we get the following expression: 

2=1 



r*~ 3 g(z) f doj 



S I 2C ? ; 



(l-z) rf " 4 R 2 



2vr 



(D.59) 



2=0 
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Considering the boundary condition in Equation D.36 we define the normalized a U)V by: 

au>.y = R 2 Aj_ (a;,v)a w , v , 

so that: 



(D.60) 
(D.61) 



Inserting this definition into the expression for the action, we get: 



where: 



J v 

ri 3 ~/ \ 



2=1 



2 = 



T±(u,v,z) = 2C V 



ral J Jz)d z a llJ ,v(z) 



(D.62) 



(D.63) 



' R d-2 (l_ z )d-4 

The components of the retarded correlators corresponding to the transverse mode are then: 



Gi L ( W ,v) = -ZF L ( w ,v,z)|_ 1 

= -4a 



i? d - 2 (1 - Z 



d-4 



2->l 



-4C V 



rl 3 1 cW v (z) 



(D.64) 



Z->1 



E Derivation of Exact Solutions for the K = — 1 Case 



In the K = — 1 case, both the massless scalar and the transverse vector equations can be 
written in the form: 



(1 - z) a d z 



9{z) 



or 



where 



g(z)d 2 z i; + 



(l-z)« 
ag{z 



d z ip 





[ A 2 2 1 


+ 


~( \ q 




_9{z) 



g'(z) + 



1 - z 



A 2 2 
9{z) 



ifj = 



a 



with the boundary conditions: 



In the case K = —1 we have 



d — 2 for scalar, 
d — 4 for transverse vector 



S(*) = l-(l-*)2. 
In order to solve the equation, first make the transformation: 

ffl E^) = l-(l- z f = z(2-z) , 



(E.l) 
(E.2) 

(E.3) 

(E.4) 
(E.5) 

(E.6) 



-45- 



so that 



d z = -^-d w = g'{z)d w = 2(1 - z)d w = 2-y/l - wd u 
az 



The equation then becomes: 

w(l - w)d 2 ip + 
with the boundary conditions: 



1 + -(a - 3)w 



A 2 _ 1 

Aw 4 ' 



Next, define 



■0 = 



(E.7) 
(E.8) 

(E.9) 
(E.10) 



where 7 = — i^, and get an equation for </>: 



u>(l - w)d 2 J 4> + 



(l + 2 7 )+ ( -(« -3) -2 7 ) to 



0„ 



+ 



7 - 7 2 + ^ 7 (a - 3) - -q 2 



with the boundary conditions: 



lu>=0 



1 



lto=l 



. 



0, (E.ll) 



(E.12) 



Equation E.12 is known as the hypergeometric equation. The solution to this equation 
that also satisfies the boundary condition at w = is 2-^1 (a, b; c; w), where: 



c = 1 — i\ 



The longitudinal vector equation can be written in the form: 



or 



where 



g(z)d 2 ^ + 



~g{z)(l-z) a d, 
ag(z 



1 



(i- z y 



+ 



500 



ip = 



1 - z 



: , 



ag{z)\ A 2 _ 2 



a = d — 4 , 



with the boundary conditions: 



^■\ z=0 ~z-c {1-zTd, 



1 



Defining again: 



(1 - z y 

w = g(z) = l-(l-z) 2 , 



. 



2=1 



(E.13) 
(E.14) 

(E.15) 

(E.16) 
(E.17) 

(E.18) 
(E.19) 
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the equation becomes: 



w(l - w)dlifj + 



Next, define 



1 + — (a — 3)u> 



d w ip 



+ 



A 2 1 2 1 

-. -A + -a + 

Aw 4 2 



aw 



4(1 - to) 



ij = w 7 (i - w) s <i> , 

where 7 = — and S = h, and get an equation for <f>: 



if> = . (E.20) 
(E.21) 



1 



(1 + 2 7 ) + f -(a - 3) - 2 7 - 25j u; 
7 " 7 2 + ^7(« " 3) - J? 2 + - 5(5 - 1) + Q(o - .">) - 2: j 5 



u>(l - w)d^(f> + 
+ 

As for the boundary conditions, at w = we get: 

C=o = , 

while at w = 1 we have 18 
1 



. (E.22) 



(E.23) 



(1 - z) a d z 



:1> 



z=l 



(1-Z) 

= 2u> 7 (l - w)^d. 



2(1 -w^du 



w • 



(1-w) 



a-1 
2 



W = l 



1 



(1 — W) 2 



+ 27(1 — w)w 



1 



(1 — W) 2 



~ 2(1 -to) 2 d u 

The solution is once again 2-^1(0, b; c; w), where: 
c = 1 — iX 



a,b=- 1 A (a-3)- l ^±y\(a + ir-qi 



w=l 



w=l 



. (E.24) 

(E.25) 
(E.26) 



F Derivation of Vector Correlators for the K — —\ Case 

For the longitudinal vector mode, we continue from the transformation defined in Ap- 
pendix E. The expression in Equation 3.7 may be written in terms of <f> and w: 



G£( W ,s) = 4CV 



r d-3 


1 q 2 s 


(1 - w)2 WiS 


R d - 2 (1 


— w) a 2 , . o+i 
1 (1 -w) 2 d w 




+ 7^^,s 




ct-l (pto.s 
.(l-to)-a- 



• (F.l) 



Ul— >1 



18 Since the exponents of the equation at tu = 1 are and 9l y^, and since d > 4, near to = 1 the solution 
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As explained in Appendix E, the solution to the EOM with an incoming- wave bound- 
ary condition at the horizon is c/> UtS = 2 F\ (o, b; c; w), where a, b and c are given by Equa- 
tions E.25 and E.26. Define: 

a — 1 d — 5 „. 

A = c-a-b = = . F.2 

2 2 y J 

There are now several possible cases: 

1. d and is even (A is non-integer). In this case, the connection formula in Equation 3.19 
holds, so that: 

4>u),s = 2F1 (a, b; c;w) = A 2 F\ (a, b; -A + 1; 1 - w) 

+ B{1 - w) A 2 Fi (c - a, c - b; A + 1; 1 - w) (F.3) 



'1 - 10) "a 1 a. 



(1 — w) 2 



(c — a)(c — b) 



Fi(o, b; c + 1; w) 



A 2 Fi(a,b;-A;l-w) 



B(l-w) A+l 2 F 1 (c-a + l,c-b+l;A + 2;l-w) , (F.4) 



where: 



,4 
B 
A 
B 



r(c)r(A) 



T(c-a)T(c-b) 

r(c)r(-A) 
r(a)r(6) 

(c-a)(c-6) r( c +i)r(A + i) 

c r(c-a + l)r(c-6+l) 

(c - o)(c - b) T(c + i)r(-A - 1) 



AA 



c r(a)r(6) 

We proceed by calculating the limit: 



Xu 



w=l—e 



2AA 
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(F.5) 
(F.6) 
(F.7) 
(F.8) 



(F.9) 



where contact terms have been dropped. Putting this into Equation F.l we get: 



4a r d + ~ 3 2 r(— A) T(o + A) T(b + A) 



d-5R d -z qs r(A) r(o) r(6) ' 



(F.10) 



The other components of the longitudinal mode correlator can be calculated from 
Equations 3.8 and 3.9. 

2. d > 5 and is odd (A > is an integer). In this case, the connection formula in 
Equation 3.24 holds, so that: 



iF x (a, b; c; w)\ w ^ = A + . . . + B(l - w) A + . . 



(F.ll) 
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(i -w) "^ 1 d u 



1 



(1 — W) 2 



w— >1 



(c — o)(c — 6) 



w— >1 



A + ... + B(l-w) A+1 + ... , (F.12) 



where: 
B 



r(c)r(A) 



r(c-a)r(c-6) 

r(c) (-i) A+1 
r(o)r(6) (A)! 



0(o + A) + V(6 + A)] 



(c-a)(c-b) r( c + i)r(A + i) 

c r(c-a + l)r(c-6+l) 



Ai4 



5 



(c-o)(c-6) r(c + l) (-1) 



A+2 



(a + A + l) + V(^ + A + l)] 



c r(a)r(6) (A + l)! 
Proceeding in the same manner as in Equation F.9 above, we have: 



Xu 



w=l—e 



2AA 



so that: 



A+l 



Gg(o;,s) = 2C 1 ,-fc< E 



(a) A (b) A [il>(a + A) + il>(b + A)] . 



(F.13) 
(F.14) 
(F.15) 
(F.16) 

(F.17) 

(F.18) 



The other components of the longitudinal mode correlator can be calculated from 
Equations 3.8 and 3.9. 



3. d = 5 (A = 0). In this case the following connection formula holds: 

<f>u,a = 2Fi(a,b;c;w 

r(c 



r(a)r(6)^ (n!)^ ^ W > Wj 



r( a )r(6) ^ [VKa + n) + ^ + n) " 2VKn + 1)] (1 " wr (F - 19) 

(while for 2 -Fi(a ; c + l;u/) the formula in Equation 3.24 still holds), so that: 

&>,-L_n = 2^(0,650; = Aln(l -«;) + ... + £ + ... (F.20) 



(1 -w)"^^ 



(1-ti;)- 



-1 fw,s 



to— >1 



(c — o)(c — 6) 



2 Fi(a, 6; c + 1; w) 



w— >1 

A + . . . + B(l - w) + . . . 



, (F.21) 
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where: 



r(c) 
r(a)r(6) 

r(a)r(6) w 

(c — a)(c — 6) 



.4 
B 
A 

B 

c 

Calculating the limit, we get: 



» + </>(&)] 

r( c + i) 



c r(c-o + i)r(c-6 + i) 

(c- a)(c- 6) r(c+l) 



-A 



r(o)r(b) 



(o + l) + ^(&+l)] 



(F.22) 
(F.23) 
(F.24) 
(F.25) 



qjB 
2 A 
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(F.26) 



w=l — e 

where contact terms have been dropped. Putting this into Equation F.l we get: 

„2 



Gf t {uj,s) = -2C v ^ql[^(a)+m] ■ 



(F.27) 



This result also matches the general result for d > 5 in Equation 3.31, so that one 
may extend that result to d = 5 as well. The other components of the longitudinal 
mode correlator can be calculated from Equations 3.8 and 3.9. 

As for the transverse vector mode, the calculations are very similar to the scalar case. 
Writing the expression in Equation 3.11 in terms of (j) and w we get: 



Gf ± (uj,v) = -8C V 



d-3 



1 



On 



R d~2 



(1 — W) 2 



(F.28) 



The solution to the EOM with an incoming-wave boundary condition at the horizon is 
again WjS = %F\ (a,b;c;w), where a, b and c are given by Equations E.13 and E.14. The 
rest of the calculations are identical to the scalar case, except with: 



A 



a + 1 d — 3 



(F.29) 



so that we get for even d: 



G I [ ± (u;,v) = 4(d-3)C v 



„d-3 



r(-A)r(q + A)r(6 + A) 

R d ~ 2 r(A) r(a) r(6) 



(F.30) 



and for odd d: 



Gf ± (oj,v) = 8C v 



r d-i (_!)A+1 

bJ 1 ^ r 2 (A) 



(a) A (b) A [^(a + A) + ^(b + A)} 



(F.31) 
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G Effective Potentials and Their Asymptotics 



In order to obtain an asymptotic expression for the QNM frequencies and corresponding 
correlation functions, the Equations 2.9, 2.11 and 2.15 must first be put into the Schrodinger 
form: 

9|0+[A 2 -F(S)]0 = O. (G.l) 
Both the scalar equation and the "transverse" vector equation are of the form 



(1 - z) a d z 



d z tp 





r \ 2 i 

A 2 


+ 


~( \ « 




_9{z) 



ip = 



where 



d — 2 for scalar, 



a 



(G.2) 



(G.3) 



I d — 4 for transverse vector 
In these cases we can bring the equation to the form in Equation G.l by defining 

iP(z) = (l-z)1cf>(z) . (G.4) 

After substituting tfj in Equation G.2 and some algebra we get the equation 

~g(z)d z [g{z)d z <j>] + [A 2 - Vi(z)] = , (G.5) 

where 

' '/(-> 



ag(z)d z g(z) a fa \ g 2 (z) 2 

2 1-z + 2\2 +l )^zy +q9{z) - (G - 6) 



The boundary conditions are accordingly 



I 2 = 



_ iA 

~ z c 



. 



z=l 



Next we define a new coordinate as follows: 

1 

z 



dz , 



(G.7) 



(Gi 



so that g{z)d z = d z . Changing to this new coordinate we obtain the form of Equation G.l 
with V\ (z(z)) as the effective potential. 
The "longitudinal" vector equation is: 



Z 



g(z)(l-z) a d z 



1 



where 



(1-z) 
a = d- 4 



+ 



A 2 



500 



V? = o 



(G.9) 
(G.10) 
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Define again 

#s) = (l-z)M*) • 
Substituting into Equation G.9 we again get 

g(z)8 z [~9(z)d z 4>} + [A 2 - V 2 (z)} = 0, 

only with the effective potential 

v 2 (z) = - ((i - z^m) 

ag{z)d z g{z) a fa 



(G.ll) 
(G.12) 



+ - I 1 

2 l-z 2 V2 



9 2 (z) 



The boundary conditions are accordingly 



I 2 = 



z~c (i- z )*d z (1-z)- 







2=1 



g^Cs) • (G.13) 



(G.14) 



Transforming again to z we obtain the form of Equation G.l with V% ( z (z)) as the 
effective potential. 

Looking at the effective potentials, we can see their behaviour near the horizon, near 
the boundary and at z — > oo (or z = Zq). 

At the horizon (z — )■ 0, z — >■ —00) we have: 



Vl(z)->0 
V 2 (z)^0 

At the boundary [z — >• 1, z — > 0) we have: 

, . a /a \ 1 a /a 

V l{ z)*-(- + l 



V 2 (z) 



a /a 



2 V2 



(1-z) 2 2V2 
1 a (a 



(1-z) 2 2V2 



+ 1 
1 



1 _ j! - 1 

z 2 4z 2 

1 = J? ~ 1 

z 1 4z 2 



(G.15) 
(G.16) 



(G.17) 
(G.18) 



so that, generally 

V(z) 



jI - 1 

45 2 



d — 1 for scalar, 

d — 3 for transverse vector, 

\d — 5 j for longitudinal vector 



(G.19) 



At the limit z — > 00 (z — > zq) we have: 



a 



Vl (z)^^(l + K) 2 

V 2 (z)~^{l + Kf 
so that, generally 

Viz] 



a 



d + 2 



z 



2d-i 



1 a 



a 



+ d-2 



z 



2d~A 



(d-2) 2 2 

1 a 

(d-2) 2 2 



o 



d + 2 



1 



o- 



+ d- 2 



1 



(G.20) 
(G.21) 



1 



4(5 - zqY 



2 

d-2 



d-2 



for scalar, 

for transverse vector, 
for longitudinal vector 



(G.22) 
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H Asymptotic Correlators for the Vector Case 



Proceeding from the definitions in Subsection 2.4, the expressions in Equation 3.7 and 3.11 
can be written in terms of <b and z: 



„d-3 



-1 



R d-2 zd-4 



„d— 3 



■l) d 5z^,v 



2-5-0 



R d-2 



(H.l) 



(H.2) 



We again continue with the method applied in Subsection 2.4 and Subsubsection 3.4.1. 
In this case, the solution around z — > zq (z — > oo) is given by Equation 2.37 in direction (1) 
and Equation 2.39 in direction (2). From the boundary condition at the horizon (z = 0, 
z — > — oo) we have: 

<f> « e~ iXi . (H.3) 
We again distinguish between two possible cases: 

1. d is even. In this case the solution around z — > is approximately: 

cf) « B+P+(z) + = B + V2ir\zJ 11 (\z) + V2vrA5J_ a (Xz) , (H.4) 

2 2 

where: 



Ji 



d — 3 for transverse vector, 
|d — 5 1 for longitudinal vector 



(H.5) 



(so that y > is non- integer) . Replacing the Bessel functions with their asymptotic 
form, we have: 



-i/3_ 



B+e if)+ + 5_e 



-iAz 



(H-6) 



where /3± = f (1 ± ji). We now match the solutions on lines (1) and (2). We have 
for line (1): 

A+e i(-Xz -a+) + A _ e i(-\z -a.) = q _ ^ 



For the section (2) we have: 



A i(-\z +3a+) _|_ ^_ e i(-A5 +3a_) _ n 



+ S_e" 



(H.8) 
(H.9) 



Equations H.7, H.8 and H.9 form a linear system of equations. Solving for B±, A. 



we can get an expression for 



B- 



ill 



2 cos (f j 



2 Jooj 



e 2ie - + 2cos(§.7 00 ) ' 



(H.10) 



where 6L = Az + f (ji - 2). 
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Next we evaluate the correlators. Developing (ft around z = 0, we again have: 



Az- A+ ^ + ... + Bz A+ ^ + ... 



with: 



.4 



D 



2- A r(-A + i 



^x 1 



-B 4 



(H.ll) 

(H.12) 
(H.13) 



ji 



(where A _ , 
we have: 



2 A r(A + 1) 

For the longitudinal vector mode (for which j\ = d—5 = a—1 = 2A), 



ji A ( 



(H.14) 



Putting this into Equation H.l, and using the fact that for the longitudinal vector 
mode joc - 

Subsubsection 3.4.1), we obtain the correlator: 



2 — ^2 an d the symmetry properties of the retarded correlator (as in 



G*(w,s) 



d- 5 



r(A) 

c 



2 /; - • 2 cos | / | 



2i6»_ 



+ 2 cos 



cZ-2 



3 2i6L 



2 cos l ^2 I <' 



■lie 



2 cos 



d-2 



(H.15) 



where: 

MA) ^ — 6T (-A*) = Xzq - |(d-7) . (H.16) 

The other components of the longitudinal mode correlator can be calculated from 
Equations 3.8 and 3.9. 

For the transverse vector mode (for which ji = d — 3 = a + 1 = 2A), we have: 



(H.17) 



Putting this into Equation H.2, and using the fact that for the transverse vector 
mode joo = an d the symmetry properties of the retarded correlator , we obtain 
the result: 



4(d - 3)C v R d 



t r(-A) /iw 



- 



r(A) V 2 



d-3 



,2i6L 



2 COS I -g±2 I 1 ' 



2i6L 



2 COS 



d-2 



3 2i6L 



+ 2 COS ( ^ I (' 



2i8- 



+ 2 COS 



<2-2 



(H.18) 



where: 



7T 



e-(\) = -e*_(-\*) = \z* --(d-5) . 



(H.19) 
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2. d is odd (so that y > is an integer). In this case the solution around z — > is 
approximately: 



B + P + (z) + B-P-{z) = B + V2ir\zJn (Xz) + V2irXzY iL (Ai 



(H.20) 



Replacing the Bessel functions with their asymptotic form, we have: 



(B+-iB-)e 



-ip+ 



e iXz + 



(fl+ + *B_)e 



-iXz 



(H.21) 



where (3 + = f (1 + ji). We now match the solutions on lines (1) and (2). We have 
for line (1): 

A i(-Xz -a+) + A _ e i(-\z -a.) = q _ 



For the section (2) we have: 



A+e i(~Xz +3a + ) + A _ e i(-Xz +3a^) = ^ + _ iB _^-iP+ 
A i(Xz +a + ) + A _ e i(Xz +a.) = , B + iB _yP+ . 



(H.22) 

(H.23) 
(H.24) 



Equations H.22, H.23 and H.24 form a linear system of equations. Solving for B±, 
A- we get an expression for 



id 



El 

B_ 



,2i0 4 



2i- 



2 COS (f Joo) 



(H.25) 



where 9 + = Xzo — f (ji + 2). 

Next we evaluate the correlators. We first assume that j\ > (this is satisfied for 
the transverse mode with d > 5, or the longitudinal mode with d > 5). Developing <j) 
around z = 0, we again have: 



Az~ A+1 2 +... + Bz A+1 i +... 



with: 



A = - 



2ir\- A+1 2 



tt2- a 



(A - 1)! B. 



(H.26) 
(H.27) 



2^A A+ I , , 

B = . /A n B 

vr2 A r A + l) V 2 



A 

27rA z> 



vr2 A (A)! 



;i) + V(a + i)] s_ + 



2^A" 



2 A r(A + 1; 



5+ (H.28) 



(where A = y ). For the longitudinal vector mode (for which j\ = d—5 = a—1 = 2A), 
we again have: 



qjB ( 
ji A 1 



(H.29) 
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Putting this into Equation H.l, and using the fact that for the longitudinal vector 
mode j'oo = 2 — an d the symmetry properties of the correlator, we obtain the 
correlator: 



e 2ie 

2iri 



4/ T FT\2 



("I) 



A+l 



r 2 (A + i) V 2 



d-5 



e 2i8+ + 2 cos 



d— 2 



2vri- 



e 2ie+ + 2 C o S 



d-2 



2 In 



iX 



where: 



7T 



.( A ) = -^(-A*) = A^ + -(d-3) . 



(H.30) 



(H.31) 



The other components of the longitudinal mode correlator can be calculated from 
Equations 3.8 and 3.9. 

Similarly, for the transverse vector mode (for which ji = d — 3 = a + 1 = 2A), we 
get: 




where: 



7T 



_( X ) = -0*_(-X*) = Xz* + -(d-l) . 



(H.32) 



(H.33) 



Next we turn to the case of ji = 0, which is satisfied for the longitudinal mode with 
d = 5. In this case, Equation 3.67 is still valid with A = Q = 0, so that: 



Az2 ln(-z) + . . . + Bzi + ... , 



with: 



A= 2V2 I X ls _ 



7T 



B 



2V2ttA2 



7T 



In 



5_ + V2vrA2 J B H 



This case is restricted for the longitudinal case, for which we have: 



Xu 



z*d s 



22 



eln(-e) + ... + -e + ... 



(H.34) 

(H.35) 
(H.36) 

g s 2 |e, (H.37) 



where contact terms were dropped in the last equality. Putting ^ into this expression 
we get: 



Xu 



B_ V 2 



e . 



(H.38) 



This result is consistent with the corresponding expression for d > 5, so that we may 
extend the rest of the results for d > 5 to d = 5 as well. 
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I Numerical Methods 



Here we outline the numerical methods used for the exact calculation of the QNM spec- 
trum and the retarded correlation functions associated with the massless scalar and vector 
perturbation modes. 

1.1 QNM spectrum calculation 

For each of the QNM Equations 2.9, 2.11 and 2.15 and their corresponding boundary 
conditions, given the values of the parameters d, K and q, one may use the Frobenius 
method to calculate the exact values of A for which the eqautions have non-trivial solutions. 
The stages of the calculation are as follows: 



1. We write the QNM equation in the form: 

Pl(z)3%i/> + P 2 {z)d z ^ + P 3 (z)rP 







(LI) 



where Pi(z), P2(z) and P%{z) are polynomials in z that depend on the values of d, 
K, q and A. We write the boundary condition at z = 1 (the AdS boundary) in a 
similar form: 

QiW + Q2WU = o, (1.2) 

where Qi(z) and Q%(z) are polynomials in z 19 . For example, for the scalar QNM 
perturbation we have: 



Pi(z 

P2{Z 

Pz(z 
Qi(z 
Q 2 (z 



~9 2 (z)(l-z) 

g(z) [g'(z)(l-z) + (d-2)g(z)] 

(1 - z) [A 2 - q 2 g(z)} 



1 . 



(1.3) 
(1.4) 
(1.5) 
(1.6) 
(1.7) 



2. We calculate the exponents of the QNM equation at z = and z = 1 via its indicial 
equation. At z = we have: 

^=o ~ > (!-8) 

where C = g'(0). In accordance with the incoming- wave boundary condition at the 
horizon, we always choose the exponent 1 = —%■ At z = 1 we get the real exponents 
(3i and fa, and we define: 

(3 = mm(fa,fa). (1.9) 

3. We make the transformation: 



V> = z 7 (i - zf4> , 

This boundary condition should be interpreted at the first order of 1 



(1.10) 

— z that doesn't trivially vanish. 



so that (j) is regular and finite at z = 0. Re- writing the equation in terms of 4> we get: 

Pi{z)d 2 z <p + P 2 {z)d z (t> + P z {z)<t> = , (HI) 
with the boundary condition at the horizon: 



Qi{z)d z 4> + Q 2 {z) 







2->l 



(1.12) 



where Pi(z), P2(z), Pz{z)-, Qi(z) and Q 2 (z) are polynomials in z. We assume that 
Pi, P 2 and P3 are given by: 



ni 

P l {z) = Y J c\z k 

k=0 

«2 

P 2 {z) = Y J c\z k 

k=0 
n 3 

P 3 (z) = J24z k 

k=0 



(1.13) 
(1.14) 
(1.15) 



4. In order to increase the numerical stability of the calculation and decrease the time 
required for calculation we minimize the degree of the polynomials in the equation 
by dividing P±, P 2 and P3 by their greatest common divisor (as polynomials in z). 

5. We solve the equation using the Frobenius method: We develop as a power series 
around z = up to the iV-th degree: 



N+l 
k=l 



-\Z 



k-1 



(1.16) 



We then put the series into the Equation 1. 11 and obtain an equation for each power 
of z up to the (iV — l)-th degree. The homogeneous system of equations reads: 



where 



Sa = , 



a 1 



(1.17) 



and 



J 'j 



\a>N j 



cl i+2 U - i)(i - 2) + c? -+i(j - 1) + c: 



-1-3 



(1.18) 



(1.19) 



(for 1 < i < N, 1 < i < N + 1 and -2 < i - j < max(m - 2,n 2 - l,n 3 )). The 
boundary condition at z = 1 gives another equation: In the case of (5\ = (3 2 , the series 
in Equation 1.16 doesn't converge at z = 1 for N —> 00, since one of the solutions 
to Equation 1. 11 goes like log(l — z) near z = 1, while the other is finite at z = 1. 
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In order to choose the finite solution (as dictated by the boundary condition in this 
case), we choose the following condition that is necessary for convergence: 



CI JV = . 

In the case of fa ^ fa, we have the boundary condition: 

hd^ + b 2 ^ = o , 

where: 



Qi(z) 



(l - z y 



b 2 



Q 2 (z) 



z-+l 



l-z) 



(1.20) 



(1.21) 



(1.22) 



2-5-1 



and s is the minimal exponent such that b\ ^ or b 2 7^ 0. In terms of the series 
expansion this gives the equation: 



where 



N+l 

3=1 



Sn+i,] = h(j - 1) + b 2 



(1.23) 



(1.24) 



We end up with a linear system of N + 1 equations, where the coefficients are poly- 
nomials in A: 



M 



S lJ i\,d,K,q)=Y,Sl l j ) (d,K,q)X 



(1.25) 



1=0 



Finally we find all the values of A for which a non-trivial solution exists to the system 
of equations, by solving the generalized eigenvalue problem: 



det(S) = det ^ Sj^AM = 



(1.26) 



6. We filter out all of the values of A that don't converge as we take N — > 00, by 
performing the above calculation for the N-th degree and for the (N + l)-th degree, 
and including only values of A that satisfy the condition: 



A (,JV+1) _ X (N) 



< e 



(1.27) 



where e is some convergence threshold. 
1.2 Correlators calculation 

For each type of perturbation discussed in this chapter, and given the values of the pa- 
rameters d, K, q and A, we outline a numerical method based on the Frobenius solution 
to calculate the exact value of the retarded correlation function of the dual gauge theory 
operators. The calculation of the correlators require the evaluation of the expressions in 
Equations 3.5, 3.7 and 3.11. The stages of the calculation are as follows: 
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1. As in Subsection 1.1, we write the relevant QNM equation in the form of Equation 1.1. 

2. We again calculate the exponents of the QNM equation at z = and z = 1 via the 
indicial equation. At the horizon we choose the the exponent 1 = —% corresponding 
to the incoming- wave condition (and the retarded correlator in the gauge theory). 
At z = 1 we get the real exponents f3\ and fai and we define: 

p = minOSi, ft) (3' = max(/3 1 ,/3 2 ) - min(^,/3 2 ) . (1.28) 

3. We again perform the transformation in Equation 1. 10, and get the Equation 1. 11. 

4. We divide the polynomials P±, P 2 and P3 by their greatest common divisor. 

5. We solve the equation around z = using the Frobenius method: We develop <f> as a 
power series around z = up to the N-th degree: 

N+l 
k=l 

We put the expansion into the equation 1. 11 and obtain an equation for each power 
of z up to the (N — l)-th degree. We end up with the system of equations: 

Sa = , (1.30) 

where S is given by Equation 1.19. To get a single solution a normalization equation 
needs to be added, and we choose: 

a?> = 1 ■ (1.31) 

Solving this system of equations for specific values of d, K, q and A, we obtain the 
solution (j>\ (which is regular at z = 0). 

6. We find the two independent solutions of the equation, (j>2 and $3 around z = 1 using 
the Frobenius method: First, we assume that: 

02 1 

We define: 

02 

and again re- write the equation in terms of 4>2- We develop 02 as a power series 
around z = 1 up to the iV-th degree: 

N+l 

2 (z) = ^4 2 2 1 (z-l) fc - 1 . (1.34) 

k=l 

Putting the expansion into the equation, we again end up with a system of N linear 

(2) 

equations for the coefficients a k (similar to Equations 1.19, but here the parameters 



(i- z y 



(l-z) 



Zf'02 



(1.32) 



(1.33) 
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ci are the coefficients of the polynomials in the QNM equation after the transfor- 
mation given in Equation 1.33, developed around z = 1). Adding the normalization 
equation: 

4 2) = 1 , (1.35) 

we solve the system of equations to get the solution 2 and then transform back to 
4>2- The other solution around z = 1 takes the form: 

3 = 03 + cln(l-z)02 , (1-36) 

where: 

03 « 1 • (1-37) 

2 = 1 



We develop .03 as a power series around z = 1 up to the iV-th degree: 

AT+1 

ta = E&-i) w > (i-38) 



and put 03 into the equation. We end up with a system of N linear equations for 

(3) 

the coefficients a) and the coefficient c. Two equations must now be added: A 
normalization equation, and another equation to fix the extra degree of freedom of 
adding to 03 a function that is proportional to 02- In the case where /3' > 0, we 
choose the following two equations: 

4 3) = 1 (1.39) 
af = . (1.40) 

If (3' = 0, we choose: 

c = 1 (1.41) 

4 3) = . (1.42) 

Solving this system of N + 2 equations we get the solution 03. 

7. We extract the connection coefficients A, B linearly relating the solution 0i to the 
solutions 02, 03: 

<t> x {z) = Afa(z) + B<fo(z) . (1.43) 

This can be accomplished by choosing a value < zq < l 20 , and setting z = z$ 
both in Equation 1.43 and its derivative with respect to zq, thereby obtaining 2 linear 
equations for A, B: 

Mzo) = AMzo) + BM*o) (1-44) 
0i(*o) = A<p' 3 (z ) + B<p' 2 (z ) . (1.45) 

Solving these equations, we get the values of A, B, and may then calculate the rela- 
tion ^ that enters into the expressions for the retarded correlators (as explained in 
Subsection 3.2). 



3 zo should be chosen so that all three series solutions are well converged in the vicinity of zo. 
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8. For the sake of better visual representation of the results, we define and calculate 
the "normalized" retarded correlators /(A) in the following way: In all cases we may 
write the retarded correlator in the form: 

G R (X) = N(X) [/(A) + M(A)] , (1.46) 

where |/(A)| = 0(1) for |A| — > oo, -/V(A) is proportional to X s (s being the appropriate 
exponent for perturbation type) and M(A) is proportional to ln(aA). For example, 
for the scalar perturbation mode (see Subsubsection 3.4.1): If d is even- 

iV (A) = -2(d-l)C s ^ r ^^y 1 M(A) = 0. (1.47) 

If d is odd- 

««= 8 »4w(?)" "w—Mt)- (l48) 

Using the numerically calculated value of ■? and the values of N and M appropriate 
for each perturbation type we may extract a numerical value for I(X). 
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